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FOREWORD 


This final technical report describes the work accomplished during the 
period from 1 September 1988 through 31 October 1989 by the University of 
Dayton Research Institute (UDRI) under Grant No. 60NANBSDO557 for the National 
Institute of Standards and Technology. Dr. Vytenis Babrauskas of the Fire 
Measurement and Research Division, Center of Fire Research was the technical 


monitor. 


The efforts reported here were performed by the principal investigator, 
Mark A. Dietenberger, a full-time member of the Structural Integrity Division 
of UDRI. Dr. A. M. Rajendran, head of the Analysis and Computation group 
provided project direction and Dr. Joseph P. Gallagher, head of the Structural 
Integrity Division served as the program manager. Ms. Joyce Smith provided 
assistance with programming expertise in the area of documentation and 
operating systems. Ms. Joanda D’Antuono was responsible for typing and 


assembling the quarterly progress reports and this final technical report. 
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TECHNICAL REFERENCE AND USER'S GUIDE FOR FAST/FFM VERSION 3 


Mark A. Dietenberger 


University of Dayton Research Institute 


Abstract 


The FAST/FFM computer model provides dynamic, quasi-three-dimensional 
predictions of furniture fire growth and burnout in a room as well as of the 
spread of both nontoxic and toxic gases, smoke, and fire to other rooms. This 
reference guide provides a detailed description of the furniture fire model 
(FFM) equations, the program structure of FFM and its insertion into FAST, the 
FFM data input requirements, and the output produced by version 3 of the 


model. 


1. INTRODUCTION 


1.1 Background 


The evolution of compartment fire numerical models in recent years has 
included new models of the fire source, of which the FAST/FFM model is one 
example. The compartment fire models, such as FAST [1]!, utilize the distinct 
gas zones concept as a compromise between a network model and a finite 
difference model. That is, basic equations describe the mass, momentum and 
energy transfer from zone to zone in a fire driven environment. The numerical 
solution of these equations were shown to be both practical and sufficiently 
detailed. There was a need, however, to develop better models of a self- 
consistent fire on a general mockup and of thermal radiation heat exchanges 


between all surfaces and gas zones. 


1 Numbers in brackets refer to literature references listed in Section 9 


The motivation for the furniture fire model (FFM) was two-fold. The 
first was to predict full scale fire tests using only bench scale fire data. 
The second was to provide acceptable predictions of the spread of fire, smoke 
and toxic gases for determining the hazard in various fire scenarios. Such an 
analytical tool has potential for reducing the number of costly full scale 
tests and for providing the fire protection community with improved predictive 
capability for fire hazard, particularly when evaluating new materials in a 


new environment. 


The FAST/FFM computer model provides dynamic, quasi-three-dimensional 
predictions of: (a) furniture fire growth and burnout in a room and (b) the 
spread of both nontoxic and toxic gases, smoke, and fire to other rooms. The 
FAST/FFM model merges FAST with the Furniture Fire Model as shown in Figure 
1-1. The figure also shows the input and output of the combined models. The 
three-dimensional aspect of the FAST/FFM model includes: mockups constructed 
of connected panels, several flame volumes attached to pyrolyzing polygonal 
bases and sides, a radiation heat exchange between facets of objects, flames, 
and gases, and the construction of multiple gas zones in a few rooms. The 
dynamic aspect of the model is a result of calculating: (a) flame spreading 
in any direction on a piece of furniture, (b) the scaled burning history of 
surface elements, (c) the varying heat, fuel, and combustion products release 
rates, (d) the temperature changes of surface elements, walls, and gas layers, 
and (e) the growth of the upper gas layer with combustion products. Almost 
any furniture material can eventually be simulated in the model because the 
model effectively scales the cone calorimeter data and the flame spreading 


data. All these features taken together make FAST/FFM unique and promising. 


1.2 Objectives and Organization 


The main objective of this reference guide is to describe: (a) the 
equations which constitute the furniture fire model, (b) the program structure 
to implement the equations, (c) the data which are used by the FAST/FIFM model, 
and (d) the procedure to use the model. Section 2 describes the structure of 
the FAST/FFM model. Section 3 explains the scaling of the physical processes 
such as the material heating, flame spreading, material pyrolysis, flaming 
combustion, and heat tranfers, from the bench scale data to the full scale 


application. Each physical process simulated in the furniture fire model is 
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explained with equations and numerical techniques in Section 4, the FAST/FFM 
calculation system. Section 5 documents the FAST/FFM program structure. The 
interfacing routines for merging FAST and FFM and the program modules of FFM 
are explained. Section 6 describes the files used by the model and Section 7 
explains the options of the model output. In the appendices, examples of 
input data files and the corresponding output are given. The example of model 
validation with the furniture calorimeter tests on furniture mockups is also 
presented. The FAST/FFM code is written in FORTRAN 77 and is intended to 
operate on the same hardware as the FAST code. However, it has only been 
tested on a CRAY Y-MP8/864 (State of Ohio Super Computer Center), DEC VAX 
11/780 (University of Dayton), and the computers at NIST. 


2. STRUCTURE OF THE FAST/FFM MODEL 


The structure of FAST model is described fully in Reference 1. Thus, 
this guide concentrates on describing the structure of the Furniture Fire 
Model and its interaction with the FAST model. The main objectives of the 
FAST/FFM development were to: (1) Improve the model predictions for the 
accelerative growth of the fire from smoldering/ignition to a peak value and 
then the gradual die-out as depending on a fire scenario; (2) Utilize only the 
cone calorimeter database and the flame spreading data as the materials input 
data; (3) Obtain an effective interface between FAST and FFM while making 
minimal changes to the FAST program; and (4) obtain reasonable computer 
efficiencies, maintain good solution stability, and construct modular program 


Links. Css 


The construction of modular components in FFM differed significantly 
from other reported approaches because of our attempt to realistically 
simulate an arbitrarily shaped fire on an arbitrarily shaped piece of 
furniture. That is, some researchers tend to treat a flame source as a 
modular program unit with the idea of combining the various types of fire 
sources within a compartment fire model. Our approach was to organize program 
modules according to physical processes and construct general geometrics as a 
common denominator among the modules instead of organizing program modules 
according to separate fire sources and giving a complete sequence of physical 
processes for each such module. Thus, we can explicitly consider: (a) 
merging of burning regions of flame spreading or thermal ignitions, (b) effect 
of gray body flames, gas zones, and the solid surfaces on the 3-D thermal 
radiation field, (c) coupling between surface pyrolysis and thermal radiation 
on a general 3-D mockup, and (d) calculating the mockup distribution of rising 


virgin temperature and burning material thicknesses. 


The coupling of physical processes between that of FAST versus FFM is as 
follows. FAST solves a system of ordinary differential equations (ODE) to 
provide the pressure, the upper and lower gas layer temperatures, the zone gas 
layer volume, and the chemical specie masses as a function of time. These ODE 
for each zone require as input the current values of the source terms, which 
are the thermal radiative heat rate, the surface convective heat rate, the 
fire plume mass rates, the venting mass rates, and the fire heat and mass 


release rates. Of these source terms, FFM provides current values of the 


thermal radiative heat rate, the fire plume mass rates, and the fire heat and 
mass release rates to FAST. On the other hand the furniture fire model 
requires as input the current values of pressure, temperature, volume, and 
chemical specie mass of each zone in the fire source room. These physical 
couplings between FAST and FFM require synchronizing the time integration 
procedures between the two models. The implementation of the FFM modules and 
the coupling with FAST as described above is hereby refer to as the FAST/FFM 


calculation system. 


Previously [3,4], we determined that to utilize the cone calorimeter 
data and the flame spread data, a model of a physical process must be analyzed 
in two phases. The first phase is the processing of the basic database to 
derive material parameters for flame spreading and to derive scaling variables 
corresponding to the cone calorimeter measured variables. The second phase 
uses the processed data from phase one as input to the predictive fomulae for 
the various physical processes. The program module for each physical process 
was carefully organized to obtain the optimal numerical solution as part of 
phase two. The structure of the FAST/FFM model thus contains two analysis 
phases. The first implements the scaling procedures for the bench scale data 


while the second phase is the FAST/FFM calculation systen. 


3. SCALING PROCEDURES FOR BENCH SCALE DATA (Implemented in FFMDAT) 


Figure 3-1 presents the organization of the FFM input database. The 
cone calorimeter data and the flame spread data as described in Reference 5 
are organized under FFM basic data. From this basic database are derived the 
databases as inputs to the various formulae of the physical processes in the 


FAST/FFM calculation system. 


The top derived database consists of the material parameters for 
ignition and flame spread. The basic flame spread data typically consist of 
time to ignition and of the flame spread rate on a horizontal sample under 
varying irradiances. The formulae for thermal ignition and flame spread rate 
are fitted to these basic data to derive the corresponding material 
parameters. The curve fitting routines are still separate from the FAST/FFM 
code because they are relatively complex and should be replaced by a simpler 
curve fitting routine. A new routine currently under development is based on 


a more effective model of flame spreading. 


The two intermediate derived databases identified in Figure 3-1 result 
from the scaling of quasi-stoichiometric heat release rate and the scaling of 
fuel pyrolysis rate as a function of burn time for a given material. Three 
constants for scaling the basic heat and mass release rates, as explained in 
subsection 3.2, can be determined iteratively with the plotting of the scaled 
data at different irradiances. That is, as explained later in Section 7, the 
plot of scaled heat release rate versus scaled burn time of different 
irradiance levels is examined to verify the coalesence of the data into a 
single curve. The resulting scaling constants are used by the FAST/FFM code 
to create the database of scaled heat and mass release rates as a function of 
material burn history. An algorithm has been developed to automate the 
calibration of these scaling constants; this algorithm has not been throughly 


tested and thus is not reported in further detail here. 


The bottom derived database contains the scaling of combustion products 
and of soot mass fractions as a function of burn time; it is explained more 
fully in subsection 3.3. The derived two soot parameters, the maximum soot 
absorption coefficient and the specific soot extinction areas, are presently 
treated as constants, but should eventually be obtained as a function of 


material burn history. These parameters are used in the formulae to calculate 
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the soot mass fraction and the soot extinction area in the fire plume as a 


function of burn time. 


Despite the inadequancies in the current data processing routines, the 
current FFM derived database is sufficient to obtain a validation of the 
FAST/FFM model with a selected furniture material (see Appendices A, B, and 
C). A more detailed explanation of modeling these physical processes during 


the first analysis phase follows. 


3.1 Ignition and Flame Spread Properties 


In our earlier report [2], it was determined that the fabric/foam 
cushion and the PMMA materials have the following time to ignition and flame 
spread rate properties. The fabric/foam cushion behaves as a thermally thin 
material and the PMMA behaves as a thermally thick material. Thus the surface 


temperature of a thermally thin material is evaluated by the formula, 


t q,(t-7) dr 
Ty select {} Bag estica, ot ao, 
while the surface temperature of a thermally thick material is evaluated by 


the formula, 


t q (t-1r) dr 
Tear Malte (ieee (2) 
fe) J/x kper 
The net surface heat flux, qs: due to net radiative heating and to 

convective cooling is the heat load on the material. The surface temperature 
response of the thermally thin material is mediated by the characteristic 
thermal thickness, pcé, while that of the thick material is mediated by the 
thermal inertia, kpc. The surface heat flux of a typical thermal ignition 
experiment usually consist of a constant radiant flux, a surface reradiation 
flux, and a convective cooling flux during the time to ignition. Thus the 


approximate analytic solutions to Eqs. (1) and (2) are respectively, 


To- Tym (a,,/h) (1-exp[-ht/(pcéd)}) , (3) 


Pee lena des Zn) lie exp [h7t/(kpc) ] erfe(h*t/(kpc)] 7) 


where the heat transfer coefficient is approximately, 
h= Petr - Ty) /(T =<.) eel ee. to etett (5) 
eq a eq a c mt ac 


By setting the surface temperature, To: to the ignition temperature and ~ 
the integrated time, t, to the ignition time, it is possible to curve fit the 
analytical solutions to the data of the time to ignition versus irradiance 
flux. The values for the ignition temperature and the effective thermal 
thickness of the thermally thin materials and the values for the ignition 
temperature and the thermal inertia of the thermally thick materials were 
obtained in this way (see Reference 2). More generally, a material could be 
porous, finitely thick or layered, which may imply neither thermally thin or 
thick behavior. Approximate analytical solutions were obtained recently for 
these type of materials exposed to constant or exponentially decaying heat 
fluxes. They can eventually be included in the furniture fire model with 


minor modifications. 


Many studies of flame spread often consider the deRis [6] classical 
solutions for the opposed flow flame spread over both thermally thin and thick 
materials due to the heating from the diffusion flame anchored at the 
pyrolysis front. Recently the effect of finite chemical kinetics was included 
by empirically adjusting the flame spread rate as a function of the Damkohler 
number [7-9]. However, we only have access to the ignition and the flame 
spread data corresponding to the air environment and may have to contend with 
non-ideal materials. At the other extreme, one approach to modeling flame 
spread is to define a surface element ignition "state" when the calculated 
surface temperature reaches ignition temperature. The heat transfer from the 
anchored diffusion flame to the heated virgin surface would still need to be 
calculated. This approach becomes impractical when we consider the small 
heating length (in millimeters) derived in the deRis flame spread model. That 
is, it is not practical to numerically construct surface elements smaller than 
millimeters in order to properly resolve the heating from the anchored 
diffusion flame. Thus, there must be new ways of applying the classical 
solutions of flame spread to our problem. Wichman [10] demonstrated a 
simplified flame spread model that gives similar results as deRis analysis. 


One result of Wichman's analysis is equivalent to a coordinate transformation, 


10 


Sym Mel G sy i) (6) 


where s is the surface distance from the flame front and Ve is the quasi- 
steady flame speed. This permits a consideration of a surface heat flux 

, with the 
8 


flame spread rate treated as constant over S.: Substituting the heat flux 


distribution from the flame front to the preheated distance, s 


distribution into Eq.(1) and solving for the flame spread rate of the 


thermally thin material, we obtain a good approximation (see Reference 2), 


7 ‘yee, 
1 OP hi rn a ae ae ; Ch) 
re (pcd) ite - T,) 


where the numerator components qk and S, are directly constructed from the 
surface heat flux distribution as shown in Figure 3-2. Likewise for the 


thermally thick material we obtain a good approximation, 


=a, /4s,/n 2 
k 
J/kpe (Ty, gee) 


These two equations have many possibilities. As one example, a 
correspondence with deRis results can be obtained by assuming a surface 
conductive heat flux as a decreasing exponential function of s. Then we 
convert the summation to an integral and set the ignition temperature to the 
vaporization temperature. The result for the gaseous surface conductive heat 


flux at the flame front is, 


ip = (encuavaye(he ty )/i2 (9) 


P 


and the exponential decay length (i.e., the thermal length) is, 


He = 2 k,/(0,0,,V,) (10) 


A generalization to deRis theory would require assuming that the 
diffusion flame is anchored to the critical ignition location in the fuel/air 


Li 


=o 
eo Giga gt TF 


Figure 3-2 Construction of ay and 5) from heat flux 


distribution, q = f(s), fixed ahead of the 
flame front at speed Ve 
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mixture rather than at the pyrolysis front. Thus the ignition temperature and 
the materials properties derived from the time to ignition data can be used 
also in the flame spread rate equations. One useful result for thermally thin 
materials is that the numerator in Eq. (7) is a constant regardless of the 
opposed air velocity unless the flame or the ignition temperatures are 


affected by the air environment. 


The other source for surface heat flux distribution ahead of the flame 
front is the thermal radiation from the plume of the anchored diffusion flame 
and from other burning regions. Equations (7) and (8) can consider explicitly 
the two separate thermal radiation sources. The thermal radiative heat fluxes 
from other burning regions can be treated as a constant during the preheated 
length, So: The radiative flux from the adjacent fire plume can be modeled as 
a decreasing exponential function of s. The thermal radiation analysis of the 
FAST/FFM model was designed to accommodate the above radiative heat flux 


distribution for flame spreading. 


Thus any flame spread data must be examined carefully for two factors. 
The first is to determine if the material can be treated as thermally thin or 
thick, or even neither. This is done by curve fitting the different formulaes 
to the time to ignition versus irradiance data and selecting the one giving 
the best fit. The other factor is to sort the conductive heating separately 
from the radiative heating in the numerator terms of the flame spread 
formulae. One way to do the sorting is to use the FAST/FFM code itself to 
calibrate the conductive heating contribution to the opposed flow flame 
spreading on a horizontal furniture cushion. This in fact was the approach 


necessary in the model validation with the furniture calorimeter. 


Another sorting approach to separate conductive from radiative heating 
is to use a flame spread apparatus (ie. LIFT) in which the radiative heating 
is well defined. That is, the radiative heating from the wall fire plume is 
miniscule at the flame front and any external radiant heating is applied over 
sufficiently long time for the surface to reach equilibrium temperatures. 

Thus only the conductive heat flux from the flame need be included in the 
numerators of Eqs. (7) or (8). The LIFT apparatus [11] has the external 
radiant heating designed to provide a distribution of equilibrium temperatures 
along the lateral dimensions of the test sample. The measurements of the 


flame spread rate, Veo along with the corresponding value of the equilibrium 


Ls 


temperature, Tom fea provide the data needed to calibrate the flame 
conductive heating term and to confirm the choice of Eqs. (7) or (8). Since 
many materials have a minimum equilibrium temperature above the ambient 
temperature in which to support any flame spread, one must be cautious in the 
applications of Eqs. (7) or (8) to the LIFT data. Indeed, the more advanced 
thermal models of flame spreading demonstrate that Eqs. (7) or (8) are only 
applicable at relatively high flame spread rates. 


The surface heat flux distribution for the upward flame spread was 
obtained from the data collected and correlated by Quintiere, Harkleroad, and 
Hasemi [12]. Equation (7) is then used for flame spreading on the fabric/foam 


cushions and Eq. (8) is used for the thick PMMA as an example. 


3.2 Scaling of Heat and Mass Release Rates (Implemented in PROCSS) 


The need for scaling the heat and mass release rates versus scaled time 
can be ascertained from Figures 3.3a and 3.4a, which are plots of the 
corresponding cone calorimeter data of a typical fabric/foam cushion at three 
irradiance levels. The furniture fire in general will start at low irradiance 
levels at ignition and then reach high irradiance levels at the peak burnrate 
at a later time. There is a real problem in interpolating or extrapolating 
the data in Figures 3.3a and 3.4a to the conditions representative of the 
furniture fire. A detailed model of the fuel pyrolysis and the heat release 
rate to address this problem is avoided for two practical reasons. The data 
shown seem to imply that some complex processes are not readily simulated even 
by the most advanced detailed model. Secondly, the detailed models take up 
too much computer time, particularly when several surface elements are used in 
the model geometry. The alternative approach is to formulate a scheme to 


rescale the data for use in the furniture fire model. 


The 1987 report [3] described successful scaling techniques for the heat 
and mass release rates, which are summarized below. The key concepts for 
scaling were to develop variables similar to quasi-stoichiometric heat of 
combustion and an effective heat of pyrolysis as a function of material burn 
history. The processing of the cone calorimeter data is as follows. The 
measurements of the chemical heat release flux, the fuel mass loss flux, the 
soot mass fraction, and the mass fraction of carbon monoxide as a function of 


time are used to obtain the quasi-stoichiometric heat release flux as in, 
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Figure 3-3b Correspond scaled heat release flux versus 
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Figure 3-4a original data of fuel release flux versus 
time for LO fabric/FR PU foam material 
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Figure 3-4b Correspond scaled fuel release flux versus 
scaled time for LO fabric/FR PU foam material 
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The net surface heat surface flux on the material sample is 


given by, 


Gs ~ 40 F Wye F IG ~ Ip F Io ve . 
where the convective heat flux and the radiative heat fluxes from the flame, 
the attenuated cone heater, and the burning surfaces are considered. The 
scaled heat flux constant parameter, q,* accounts for systematic errors in the 
other heat flux terms and for other unaccounted physical processes. Such 
physical processes may include glowing combustion, radiation blocking by 
stagnant pyrolytic gases, or altered pyrolysis processes at elevated 
temperatures or irradiances. The scaled heat release flux is equivalent to 
the quasi-stoichiometric heat of combustion divided by the effective heat of 
pyrolysis. The scaled mass release flux is equivalent to the inverse of the 


effective heat of pyrolysis. These are given by the equations, 


ox e 
ho./4, ils 9) 
mn = m"/q (14) 


The scaled time (or the burn history) was found to be best represented 


by an effective heat release per area as given by the following equations. 


* t e 
T = J hi, (1+ G6) de (15) 
t e e 
G=0 for J h., dts : Qep G/(1 + 4) 
or = Q. -e Oona + d) (16) 


G= (Qof Pi Q,) /(Q, se Qeef? 
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Q, - f h. dt (17) 


(18) 


G is a time stretching parameter that ensures the scaled burnout time is 
identical for the similar materials. These materials can have variations due 
to non-ideal material construction or to variability in foam thicknesses. 

Most likely these similar materials will have nearly the same burning behavior 
at the beginning and the ending part of their burn history. Thus the middle 
part of a material’s burning history will be only affected by the time 
stretching parameter as controlled by the scaling parameters c and d. 


Reference 3 give a more detailed explanations of these parameters. 


The results of using the scaling Eqs. (11) to (18) are shown in Figures 
3-3b and 3-4b, which demonstrate the approximate collapse of the data into a 
single curve after optimizing the three scaling constants. Thus the data 
presented in Figures 3-3b and 3-4b permit arbitrary variations in the net 
surface heat flux and in the burning history in order to calculate the heat 


and mass release flux of the surface elements in the furniture fire model. 


3.3 Scaling of Combustion Products and Soot (Implemented in PROCSS) 


The cone calorimeter measures currently the amount of fire products such 


as cO,,, co, H, 


time. The soot extinction area as a function of time is determined from the 


O, HC, HCl, and soot with sampling techniques as a function of 


soot absorption coefficient, the volumetric flow rate, and the fuel mass loss 
rate as a function of time. The soot absorption coefficient is measured with 
a laser attenuation technique and the volumetric flow rate is obtained with 
the temperature and air velocity sensors. Our example database is from an 
earlier version of the cone calorimeter, which only derived accurate values 
for the heat and mass release fluxes, the oxygen consumption rates, the 
product ratio of CO to CO. , and the soot extinction areas as a function of 
time. In many cases the soot extinction area data were available only in 
certain forms, such as the averaged value or the peak value. The other 


combustion products were estimated assuming stochiometric conditions within 


the fire, as the water vapor for example. With these available data the only 
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two additional data plots of significance are the c0/CO,, ratios and the soot 


extinction areas versus time as shown in Figures 3-5a and 3-6a, respectively. 


The combustion products, excluding the soot, in the current model are 
scaled the same way as the mass release flux, i.e., through equations similar 
to Eq. (14). It is therefore assumed that Eqs. (13) and (14) and those for 
the combustion products are only a function of the burn history and of the 
oxygen concentration in the air. The scaling for soot production is different 
in the sense that it is also affected by the size of the flame and by the 
chemical heat of combustion for the fire. This is given by the following 


equations from Reference 4. 


e e 
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It is assumed that the specific soot extinction area, on and the 
maximum soot absorption coefficient, Aa eee are only functions of the 
material burn history and of the oxygen concentration in the air. These two 
soot parameters are temporarily treated as constants for our example cone 
calorimeter data. Equation (21) is equivalent to Bard and Pagni’s [13] 
expression for the flaming soot absorption coefficient and was derived in 
Reference 4. The exiting soot extinction area is calculated with Eq. (20) , 
which was derived in Reference 4. Equation (19) provides the soot mass 
fraction by definition and is plotted in Figure 3-6a. The scaled soot mass 
flux, which is the product of Eq. (19) and (14), is plotted versus the 


material scaled burn time in Figure 3-6b. 


Since the production of soot within the fire also affects the resulting 
chemical heat release rate of the fire which in turn determines the flame 
size, there is a certain degree of coupling between the variables of physical 
processes that must be accounted for in the model. The calculation of soot 


quantities in different forms is very significant. In large fires, the soot 
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dominates the thermal radiation, which in turn dominates the fuel pyrolysis 
and the consequential burning rate. The smoke leaving the fire is the major 
factor in hazard analysis. That is, the hot smoke can lead to human 
incapitation, obscuration, surface damage, and room flashover. The presence 
of smoke is a common means of fire detection. Thus any model that describes 
the combustion processes in the fire source of a compartment fire model should 
include reasonable calculations of soot production. The scaling of soot 
production is therefore viewed to be as important as (1) determining the 
ignition and flame spread properties and as (2) scaling the heat and mass 
release fluxes. The scaling of other combustion products could be kept simple 
in comparison to the scaling of soot production because of their lower effect 


on hazard analysis. 
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4. FAST/FFM CALCULATION SYSTEM (Implemented in FFMGO and FFMINT) 


As noted above, there is a three fold coupling between the physical 
processes of fuel pyrolysis, soot production in the fire, and the thermal 
radiation due to the presence of the fire. It is also clear that additional, 
but different, couplings are involved in the thermal heating of virgin 
surfaces, the spreading of the fire base, and the burning history of a surface 
element. The thermal radiation field is the common denominator involved in 
these two sets of couplings between all physical processes. As a consequence, 
special attention was given to developing a thermal radiation analysis that 
could accommodate all the modeled physical processes and also be tolerant of 
the couplings between the variables. A modified form of the Hottel’s thermal 
radiation model reported in Reference 2 was successfully implemented for a 
general mockup and fire, but the coupling variables had assumed fixed values. 
As these coupling variables began to be calculated as a function of time in 
the furniture fire model, the solution accuracy, stability, and efficiency 
became a serious problem. This problem persisted in the first version of 


HEMFAST. 


The problem was resolved by analyzing the nature of couplings between 
variables describing the physical processes as shown in Figure 4-1. It was 
first noted that the virgin surface temperature response, the flame spread, 
and the burn history all changed slowly enough with time and were weakly 
coupled to each other such that they could be integrated explicitly with time 
(top part of Figure 4-1). The predictor equations for these variables were 
reformulated to be more tolerant of the time step size. The other set of 
coupling variables, i.e., the fuel pyrolysis rates, the soot production, and 
the flame size along with the thermal radiation are strongly interacting and 
involve short time processes. These variables are solved in a convergent 
iterative scheme as shown in the lower part of Figure 4-1. Finally, the FAST 
output variables tended to have a relatively small influence on the furniture 
fire, so that in a particular time step, the furniture fire growth algorithms 
are executed first so as to have an immediate impact on FAST’s predictions for 


the room environment. 
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Figure 4-1 FAST/FFM calculation system 
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4.1 Furniture Geometry Construction (Implemented in MOCKUP) 


The implementation of the models of thermal ignition and flame spread in 
the furniture fire model is based on the geometrics chosen for the model. 
Three-dimensional vector algebra is fully utilized and is even required by the 
thermal radiation model to describe the geometries for the mockup, flame 


spreading regions, flame shapes, room shapes, and the gas zones. 


A general but practical geometrical construction for the mockup is to 
simulate major mockup surfaces by trapezoidal panels which contain square-like 
surface elements. An example of a square panel with nine square elements to 
represent a mockup cushion is shown in Figure 4-2. The mockup shape is 
simulated by assembling the panels in the following way. First, a panel is 
defined in a neutral planar coordinate system. The corner points of the panel 
are defined in a counter clockwise direction around the normal vector of the 
panel using inputs of panel’s length and angle specifications. The additional 
input value of element area results in specifying the surface elements 
midpoint vectors, normal vectors, and actual areas which are the required 
geometries for the thermal radiation analysis. The user information input for 
the three-dimensional translation and rotation results in putting the panel in 
any place and in any orientation in the room. A mockup or even a number of 
mockups are constructed by repeated applications of the previous steps for the 


panel constructions. 


4.2 Prediction of Surface Temperature, Ignition, and Flame Spread 

(Implemented in TIFSST) 

The surface temperature and thermal ignition are predicted at each 
virgin surface element midpoint by the appropriate numerical integration of 
either Eqs. (1) or (2). The varying net surface heat flux is obtained from 
the furniture fire model rather than prescribed by the cone calorimeter. The 
net radiative heat flux to the element midpoint is obtained from the thermal 
radiation model and the convective heat flux is obtained from the empirical 
relationships for natural convection. The numerical integration of Eqs (1) or 
(2) can be done by a predictor-corrector scheme, but unacceptably small time 
steps may be required to maintain accurate solutions due to the surface heat 
flux also being a function of the surface temperature. Thus the predictor 


equation for Eq. (1) was devised to exponentially approach the equilibrium 
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temperature for a constant irradiance level at any time increment chosen (see 


Reference 4). The result is, 
T (ttAt) = Dee (T(t) - Seva exp[-h At/(pcé) ] (22) 


nee nay /bl eat (a(t) reDET! (te) Lay h T(t) Jn (23) 


where the approximation to Eq. (5) is given by, 
h=h+h =4ce0T (ty? +h (24) 
res c Ss c ’ 


and q(t) is the net radiative surface heat flux at current time. A 
correction to h. can be obtained by resubstitutions of den from Rqas (25) eincto 
Eq. (5). The Eq. (22) can also be used as a corrector integration when the 
equilibrium temperature, ee is a function of q.(ttAt), To (t+dt), and 
T(t+At) instead. This integration scheme will prevent a very thin material 
from obtaining unreasonable surface temperature values, thus removing one 


source of solution instability in the FAST/FFM calculation system. 


The numerical integration of Eq. (2) for the thermally thick materials 
is more difficult because one must account for spatial variation of the 
temperature within the material. Since the surface heat flux, qe? imiEqh 362) 
cannot be formulated as an explicit function of time in the furniture fire 
simulations, it is formulated as an accumulated stepping function. The 
Harvard Computer Fire Code [14] has a numerical scheme that effectively 


addresse these difficulties with the result, 


T (teat) = T. + (T(t) - T.)/(1 + 9.65 p(tat) at)’/* + 
7 ae ; 1/2 


[a,(t) + q.(teat)] + [At/(mkpc)} 3 (25) 
p(t+At) = q(trat)/fo a. (1) ar for T.<T,+ 5 and (26) 
p(t+At) = In[q,(t+dt)/q.(t) ] Ju&t for TT. Tt SE (27) 


To apply these equations to the furniture fire model, we initially set 
q,(ttAt) to q(t) so that Eq. (25) predicts Ti (t+dt). Then the net surface 
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heat flux for the next time step is calculated and Eq. (25) can be used to 
correct the predicted surface temperature. Note that the surface temperature 
increment is proportional to the square root of the time step; it is more 
likely the surface temperature values will be reasonable. In any case, the 
limit on the surface temperature value is the equilibrium temperature, Cre 
obtained from Eq. (23). Equations (22) to (27) was implemented in the 
subroutine SURTMP of the FFM code. 


When the predicted surface temperature at a surface element over a time 
step exceeds critical temperatures, the model will redefine certain parameters 
to initiate: (a) the material burn history, (b) a burning region, or (c) both 
at the surface element. If the critical temperature is the ignition 
temperature, then the following is implemented. The surface temperature is 
reset to the burning temperature, then the remaining burning time in the time 
step is calculated, the surface element is labeled as burning, and the flame 
spread region is defined with the surface element. If the critical 
temperature is the pyrolysis temperature and less than the ignition 
temperature, then the remaining "burning time" in the time step is calculated 
and a pyrolysis state is defined for the surface element. The definition of a 
pyrolysis state allows for the possibility of a material to decompose into 
gases not involved in flaming combustion at the surface element. If a surface 
element is already labeled burning or pyrolyzing and also the net surface heat 
flux, qe? becomes negative, then the burn history is frozen and a surface 
temperature decrement over the time step is calculated instead. The logic 
described above is implemented in subroutine INTIAL of the FFM code. 
Conceptually, a fire suppression and eventual re-ignition of a surface element 


can be analyzed through the effects on the net surface heat flux. 


The opposed flow flame spread region is by definition the time 
integration of the flame spread rate in the direction normal to the flame 
front and along the surface. The assisted flow flame spread region on the 
other hand spreads along the fire plume direction. An octogon with its 
vertices moving outward as shown in Figure 4-2 seems to be a practical 
representation of the flame spread region originating at a piloted ignition 
point. An elongated eight sided polygon would be practical for a line 
ignition source. Each vertex of a polygon is represented by a three- 


dimensional vector sequenced in a counter-clockwise direction around the 
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polygon normal vector. Thus the direction of the flame spread vector is a 
normalized cross product of the polygon normal vector with the vector between 
Staggered vertices. The magnitude of the flame spread vector is then the 
flame spread rate multiplied by the time step, unless it is cut short by 
intersection with the panel edge. The logic described above is implemented in 


the subroutines, FIRPNT, FIRPLY, and SPRDLN of the FFM code. 


For the prediction of the flame spread rate, the preheated temperature 
and the heat flux distribution shown in Figure 3-2 are the additional 
parameters needed in Eqs. (7) and (8). The preheated temperature is taken to 
be the predicted surface temperature at the virgin surface element midpoint 
which is the closest to the vertex of a flame spread vector. The surface 
element midpoint is then extrapolated (parallel to the flame front) to the 
flame spread vector line to define the preheated length, So: The external 
radiative heat flux determined at the surface element midpoint is also 
extrapolated as a constant heat flux as part of the heat flux distribution in 
Figure 3-2. The radiant and the convective heat fluxes between the heated 
surface and the ambient air are included in the heat flux distribution by 
treating them as a product of a linearized heat transfer coefficient and the 
temperature difference between the air and the surface. The increase of 
external radiative heat flux from the distance S, to the vertex point at s=0 
due to the adjacent fire plume is modeled as a decreasing exponential function 
of s. Lastly, the conductive heat flux distribution from the anchored 
diffusion flame was assumed constant over the distance Le. All these heat 
flux distributions are substituted into Eqs. (7) or (8) and the summations are 
converted to integrals. The Duhamel theorem and the Laplace tranform are 
applied to obtain an explicit analytical expression for the flame spread rate 


with the respective results, first for the thermally thin, 


3 T. = [a,2,+ VeeSet TPP eg npr e728) | he | poco Vel + 
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and then for the thermally thick material, 
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These two equations can be inverted algebraically to solve explicitly 
for the flame spread velocity, Ve. The conductive heat flux distribution was 
not split into two parts as was done for the external radiant heat flux 
distribution for the following reasons. Its effective length, fe, is 
typically much smaller than a surface element size, which is the same 
magnitude as S, on the average. Thus the conductive heat flux from the 
anchored diffusion flame barely affects the preheat temperature at the virgin 
surface element midpoint. Secondly, the source of this conductive heat flux 
is only from the direction of the adjacent flame front and tends to be well 
defined. Also, a minimum practical value of the flame spread rate is needed 


to initiate flame spreading from a point ignition source. 


The external radiative heat flux on the other hand tends to be long 
range as well as affected by other fire sources. Furthermore, the radiative 
heat flux tends to be miniscule at the beginning of the fire spread. The 
inclusion of the external radiative heat flux in preheating as well as in 
flame spreading should reduce the sensitivity of the flame spread rate to the 
surface element size and provide a approximate separation of the sources of 
external radiative fluxes. Note that Eq. (29) was designed to be consistent 
with the flame spread concepts associated with the LIFT apparatus, although 
there is the questionable assumption of treating the conductive heat flux from 


the flame as a constant over the length 2 Thus one might utilize the 


f° 
material thermal constants derived from conventional data processing of the 
LIFT data while understanding the inherent limitations of using Eq. (29). The 
flame spread Eqs. (28) and (29) and the associated logic were implemented in 


subroutines SPREAD and FLMSPD of the FFM code. 


The remaining features of the flame spreading procedure included in the 
model are the capability (a) to merge overlapping burning regions (implemented 
in PMERGE) and (b) to create new burning regions in the adjoining panels as a 
result of flame spreading, thermal ignitions, or even new piloted ignitions 


(implemented in SPRDLN). The three-dimensional vector algebraic analysis was 
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again applied to implement these features. This analysis ensured the flame 
spread polygons are placed properly with their corresponding panels. Once all 
such polygons are in place and all previous surface elements parameters are 
updated and redefined, the burn histories of burning elements are calculated 


next. 


4.3 Prediction of Local Burn History (Implemented in BNSTEP) 


As described in paragraph 3.2 and demonstrated in Figures 3-3 and 3-4, 
the burn history of a furniture material was sufficiently complicated that a 
practical prediction of material burning required rescaling the cone 
calorimeter data. Although Eqs. (11) to (18) have proven adequate for 
modifying the burn rate and the burn time for collapsing the data unto a 
single curve, it was not obvious how Eq. (15) can be used for the burn history 
in the furniture fire model. The problem is made more clear by taking the 
scaled heat release flux as a function of the scaled burn history from Eq. 


(13) and substituting into Eq. (15) to get, 
ok * 
Ligeia Aaa Gah (T), (LEE QCD dy yiidt (30) 


Several difficulties can be found in the use of this equation. For one 
thing the scaled burn time is an integral-non-linear function of itself as 
well as being a function of real time. The integral in Eq. (15) is evaluated 
with a trapezoidal numerical integration to produce the scaled data in Figures 


3-3b and 3-4b. A practical alternative to Eq. (30) is its rearrangement to, 
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t; T; 
The right term of this equation utilizes the scaled data presented in 
Figure 3-3b and the left term uses the net surface heat flux derived in the 
furniture fire model as a function of time. The selection between the three 
data sets in the figure (i.e., circles, triangles and x's) is determined by 
the magnitude of the net surface heat flux. Thus Eq. (31) provides a one to 
one relationship between the real burn time, t;> and the scaled burn time, Dp 


The implementation of Eq. (31) in the FAST/FFM calculation system involves a 
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trapezoidal numerical integration of the right side of Eq. (31) and an Euler 
numerical integration of the left side. Thus it was possible to keep the 


local burn history in step with the local real burn time in a simple way. 


The receding depth of the burning material is assumed approximately 
proportional to the scaled burn time. The proportionality constant is the 
virgin material thickness divided by the ending scaled burn time (or the total 
quasi-stoichiometric heat release per area of the material). If all the 
surface elements in a panel are burning, then the flame spread polygon is 
shifted to keep in planar contact with a nearest burning surface element. 

When the flame spread polygon is finally shifted by an amount equal to the 
virgin material thickness the panel no longer supports pyrolysis and 


combustion. 


The predictions of the surface temperatures, the flame spread polygons, 
and the burn histories are in essence the integrations of first-order ordinary 
differential equations which are assumed weakly coupled and non-stiff. 
Explicit integrations of these equations have not presented any problems to 
date. If the coupling or stiffness of these equations become important, the 
FAST/FFM calculation system should then consider the FAST method of solving 
the stiff system of ODE’s or other alternatives [15]. The formulae that 
describe the combustion rate, the flame shapes, the soot production, and the 
thermal radiation on the other hand are a system of coupled nonlinear 


functions. Their solution is described below. 


4.4 Local Convective Heat Flux and Fuel Pyrolysis Rate 

(Implemented in THMOUT) 

An effective solution scheme developed has all the relevant variables 
except for the radiative heat fluxes involved in the fixed point iteration 
method. The values of the relevant variables are input in an effective Gauss- 
Seidel matrix solution of the Hottel’s thermal radiation equations. The 
radiative heat fluxes thus obtained are used to modify the values of the 
relevant variables which are fed back into the thermal radiation equations. 
This process is repeated until convergence is obtained. The first step in the 
iteration, as shown in the middle of Figure 4-1, is the local evaluation of 
the convective heat flux and the fuel mass loss rate. The coupling between 


the fuel mass loss rate and the convective heat flux as well as the decoupling 


us 


between the fuel mass loss rate and the burn history is better understood by 


rearranging Eq. (14) as, 
ae i A ox r* 32 
mer Ady td de) em) CL) ; (32) 


Since the scaled burn time has been obtained from the solution of Eq. 
(31) and the scaled mass loss flux (i.e., inverse of effective heat of 
pyrolysis) hardly vary with changes in the net surface heat fluxes, the values 
of qo and qj. can be changed arbitrarily during the iterations to provide 
updated values of the fuel mass flux, m". The value for the radiative heat 
flux is assumed available from the latest solution of the Hottel’s thermal 
radiation equations. The convective heat flux is, however, a non-linear 
function of m due to the blowing effect of the fuel pyrolysis. This 
complicates the solution of Eq. (32) for the fuel mass loss rate. A very 
efficient solution is described in the following for both horizontal and 


inclined burning surfaces. 


The formulae for the convective surface heat transfer in pool fires was 
taken from Orloff and deRis [16] and in wall fire from Ahmad and Faeth [17]. 
How these formulae are used in the furniture fire model is significantly 
different than their use in References 16 or 17 or in any other references. 


Starting with the pool fires, the convective heat transfer is, 


h_(T 


Ic, pool” Dest 


Sa /SLeAP\Y Jezel). } (33) 
where 
-m"c /h (34) 
2 ae ly 
and from scaling of the fuel pyrolysis rate, 
m= = (q,+ 4, +49) / He (35) 
where h_ is the effective heat transfer coefficient in a pool fire, given in 


Reference 16 as 0.008 yey Ty 


Kelvin; Coa is the heat capacitance of air, and He is the effective heat of 


is the flame temperature, taken to be 1250° 


oe 


pyrolysis as a function of burn history. Since the convective heat flux term 
appears in Eq. (33) and in Eq. (35), the result is a nonlinear equation in the 


dimensionless mass flux y, as, 
h (Te Tey lexp Cy) eo a eee (hHe y/e,,) = 0 (36) 


Fortunately, the left hand expression is a decreasing monotonic function 
of y and is bounded. The minimum and maximum bounds of y can also be directly 
calculated. Thus the slope intercept method of solving Eq. (36) should be 
rapidly and absolutely convergent. After solving for y, Eq. (34) is used to 
obtain a value for the fuel mass flux, and Eq. (33) is used to obtain a value 


for the convective heat flux in the pool fire. 


The convective heat transfer for wall fire from Reference 17 is 


Io wall mvt ns Se 
m" = maximum(m" , m* ) (38) 
p 3(Btr)n t+ 7 
Be eo) 8 Sey ar Bon Glee) ai ceiesee sect ate 1/4 
my tk0 625 Te) ete Lancers eee besryzyerry ce mee ae 
p 3(Bt+r)no+ 7 
ae 2Oiy ei Bae 1+8) 2n(1+ 7> i 270 
me= 0.051554 (5*) (745) arr a i pesparrersyce (40) 
where 
p= coe (Tg TE /aLes, (41) 
€ = 0.5 from data correlations, (42) 
PNCL TORY 23.2% (43) 
ner lL - (rg /Ber gh yr, (44) 
Trp = 0.19 for polyurethanes , (45) 


and the local Grashof number at inclined angle, ¢, is, 
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3 2 
Gr = Les x cos(¢) / 4 Shei , (46) 
Combining Eqs. (35) and (37) we obtain a nonlinear equation for Ly as’ 
(Hp- Le) - (a,+ 45)/m" = 0 (47) 
where m" is also a function of Le via Eqs. (38) to (46). The left side 
expression of Eq. (47) is a decreasing monotonic function of L. and is 
bounded. The variable L, is bounded between O and H,. Thus the numerical 


Le f 
scheme for solving Eq. (47) is the same as that used for Eq. (36). Since we 


are able to obtain values for L- and the fuel pyrolysis rates, Eq. (37) is 
used to obtain a value for the convective heat flux on the wall. 

The determination of the convective heat flux over burning surfaces at 
different orientations is important during the early phase of the furniture 
fire growth, because convective heat from the flames then is a major factor in 
fuel pyrolysis. Since these equations are involved in a larger iteration loop 
to determine the radiative heat fluxes, it was essential to develop the most 
efficient method of solving Eqs. (36) and (47). The values for the net 
surface heat fluxes are also inputs for determining the mass and the heat 


release rates of the adjoining flame as follow. 


4.5 Chemical Heat and Mass Release Rate of Adjoining Flame 

(Implemented in THMOUT, TOFIRE, PLFIRE, WLFIRE) 

Before considering the contribution of a burning surface element to a 
flame, a set of burning regions adjoining a single flame must be determined. 
In addition, several such flames need to be determined until all burning 
regions are accounted for. A practical scheme was developed to organize the 
construction of multiple flames. The first step was to find a group of nearly 
horizontal flame spread polygons on which to form pool-like flames. If a 
nearly horizontal polygon is connected to another flame spread polygon in an 
adjoining horizontal or vertical panel, then the two flame spread polygons 
should anchor a single flame. It is possible for a set of such connected 
polygons to anchor a single flame. Such a flame is described as a merged 
pool/wall fire. However, if the derived flame height is below the height of a 


connected vertical polygon, then a separate wall fire is considered for the 
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vertical polygon and the original merged pool/wall fire is no longer anchored 
to the vertical polygon. After examining all the nearly horizontal polygons 
for the construction of merged pool/wall fires, all the remaining flame spread 


polygons will then each have a simulated wall fire. 


This organization of multiple flames construction led to two basic flame 
shapes. The first is the parallelepiped flame shape with the horizontal flame 
spread polygon as a base for the merged pool/wall fire. The other is the 
wedged-like shape defined by the Ahmad and Faeth boundary layer model [17] for 
the wall fire. The thicknesses of a wall flame for the laminar and the 


turbulent flow are respectively given by, 


3(B+r)n.+ 7 
on (1+8) £n(1+8) 41/2 £ -1/4 
Ore Cae ose T°” [yype7aciepyetx! * oe 
3(Bt+r) n+ 7 
a (1+B) £n(14B) 4/5 f -1/10 
6- 0.199571 x [5 ", iP [ayaa tx! (49) 


As a consequence, only one additional parameter, the solid flame height, is 


needed to define the overall flame shape. It is summarized below. 


Cetegen, et al, [23] formula for averaged visible pool fire height is, 


a3 (0 aie Oa 
2, /D = (50) 
P 
wk Wee i 
where 
ox ° 1/2 «5/2 
Q) > (8./10,6,7,87° 0”) (51) 


Hasemi’s formula for averaged visible wall fire height is, 
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Z {2s = (52) 


where 


Q ee (Q./Le Gael nee Z nh vis) } (53) 
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Interpolation of flame height for the merged flame is; 


Q. She 
I ony. —2 iz caaad (54) 
2 ne 6 tO ° 
P Cc Q. 
with a solid flame height: 
veh = 0.4 Z (55) 


Both flame shapes have surface elements with their positions, areas and 
normal vectors defined so that view factors for the thermal radiation model 
can be numerically calculated among the flames and the mockup surfaces. Once 
a set of burning surface elements are assigned to a flame, the fuel mass 


release rate into the flame can be calculated by a summation, 


M 5 n(T-) A (56) 
rastatn Ci rbot. "7 tka ae Bead pars beni aia 
and the combustion products mass rate is given by, 
e oX* * * 
My = y (q, + 4, + 49)5 (Ty) F(T) Ay (57) 


The chemical heat release rate for the flame is not as straight forward 
because of the need to account for the incomplete combustion due to the 


production of soot and carbon monoxide. This is modeled by the equation, 


37 


Q 24 p Gh us qr Hy qo) i [a" m (F, oe ‘ ss Hoo) litt = ey 
where E-6 is the mass fraction of carbon monoxide as a function of the scaled 
burn time of the surface element, i, and F. is the soot mass fraction of fuel 
calculated from from Eqs. (19) to (21) using mass averaged values for o. and 
: maxtor the flame. The chemical heat release flux of the fuel from a 
BUrEAce element is calculated as the quasi-stoichiometric heat release flux 
minus the heat release flux of the incomplete combustion products of soot and 
carbon monoxide. The chemical heat release rate of the flame is then the 
summation of the chemical heat release flux times the area of the burning 
surface elements. Equation (58) provides a framework for the contribution of 
each predicted physical processes to the heat release rate of the flame. The 
virgin surface heating and the flame spread ultimately determines the number 
of burning surface elements as a function of time and are a major factor in 
the structure of the thermal radiation field. The scaled burn times of the 
surface elements directly determine the scaled heat and mass release fluxes 
and the carbon monoxide mass fractions. The soot specific extinction area and 
maximum soot absorption coefficient, which are inputs for the soot mass 
fraction formula, are generally functions of the scaled burn time but are 
constants for the present version of FAST/FFM. The flame heat release rate 
itself determines the flame dimensions, which in turn affects the soot mass 
fraction (see Eq. (21)) and are also a major factor in the structure of the 
thermal radiation field. Examination of Eq. (58) also suggest the 
accelerative growth of the flame heat release rate is due to two main 
competing factors: the rate of increase of number of burning surface elements 
and the rate of increase of the surface radiative heat flux. Eventually, a 
peak value of the flame heat release rate is reached and then its rate of 


decrease is controlled by the burnouts of the surface elements. 


4.6 Thermal Radiation Analysis (Implemented in FLAMES) 


The last step in the fixed point iteration scheme is the evaluation of 
the thermal radiative heat transfer among the mockups, the room, and the gas 
zones. The thermal radiation analysis is the foundation and also the most 
complex submodel of FFM. In large fires the thermal radiation dominate nearly 


all physical processes such that it is worthwhile to have detailed geometrical 
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features and to include coupling effects between thermal radiation sources. 
Since only a few gas zones and not very many surface elements are considered, 
the Hottel’s zonal model was chosen for the thermal radiation analysis. It 
was essential to also develop a practical model, because the thermal radiation 
analysis must be done at each iteration step during each time step. The mean 
beam length approximations were then applied to the key components of the 
Hottel’s zonal model and the results were reported in Reference 2. One result 
was to simplify the evaluation of the complex Hottel’s exchange areas between 
the flame and the surfaces to only involve calculations of the configuration 
factors, surface areas, and the emissivities. An another result was the 
discovery of a very efficient scheme of solving the thermal radiation 
equations by taking advantage of energy conservation requirements. Since 
complex shapes are generated for the flames and for the mockup, the 
configuration factors were evaluated by numerical integrations. This in turn 
has made the Hottel’s exchange area calculations by far the most time 
consuming of all the routines in FAST/FFM. These calculations, however, are 
essential for accurately determining the flame spread and the burning rate on 


the mockup as a function of time. 


4.6.1 Hottel’s Zonal Thermal Radiation Analysis 


The zonal radiation analysis begins with the Hottel’s zone method 
[18] for an enclosure with isothermal, optically thin, and elemental gas 
volumes. This method was extended to isothermal, moderately optically thick, 


and finite gas volumes by considering the following. 


The surface energy balance equations: 


w, - €,e, + p;h, (59) 
iJ n co @ e 
A. hy? rhsyil(sssijwyet) pat (pes.)w (60) 
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The volume energy balance equations: 


os 16) ee et Pes (61) 
Y af af af R 


39 


n 2 
taVh-= pg )w. + ‘ 62 
ayy Be (s58,)¥; a (6,8), Se 


The net radiative power equations: 


Qo ig 7 Aga 7 AGGhy- 4) (en) 
re = an ease “ey (64) 
where 
Ww, is the outgoing surface radiative heat flux yom 
h, is the incoming surface radiative heat flux (ean) 
25 is the blackbody surface radiative heat flux (oT) Ay oe 
Sake is the surface net radiative source power (kW), 
A; is the surface element area tna)" 
és is the surface element emissivity, 
Pi is the surface element diffused reflectivity, ‘ 
Wy is the outgoing elemental volume radiative heat flux Noa tok 
hy is the incoming elemental volume radiative heat flux eee ) ; 
°, is the elemental volume blackbody radiative heat flux (gre ) (kW/m) 
oe is the finite volume radiative source power (kW), 
2. is the volume scattering coefficient, 
2 is an effective absorption Satie ak (1/m), 
Ab is an isothermal finite volume (m ) 


‘ 2 
S.S, is the exchange area between surface A. and the surface A,, (m ) 
. 2 
s, is the exchange area between the volume the and the surface A, (m), 
and 


gg is the exchange area between the volume Ms and the volume a Ee 


The SEALE Bch and g,8., are the exchange areas which embodies the 
furniture and flame geometries and the transmissivities, Tr In order to more 


conveniently evaluate these exchange areas we have the following reciprocity 


relationships, 


S.S, = S.S. (65) 
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d a 
an &,85.> 8.8 (67) 


Under the conditions of thermal equilibrium and with the Kirchkoff's 
relation e+ Pp, 7 1 for diffused emission, the energy conservation requires 
that Q- ah = 0. Substitution of this condition into Equations (63) and (64) 
and from there into Eqs. (59) and (61) result in wi ee hy and wo e, = h,. 
Substitution into Eqs. (60) and (62) then states the following energy 


conservation law: 


n Q 
AL =) 0.8.5, Fee mees (68) 
sh j=1 al eae Tal Th 
and 
n Q 
4a V =, ) s.g GREY. fee (69) 
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If the surface emissivities are less than unity and the surface 
reflectivities are finite, then Eqs.(59) to (62) can easily be solved by 
matrix techniques for the w's and h's. Equation (63) is used to obtain the 
net surface radiative power and Eq. (64) is used to obtain the net volume 
radiative power. The conservation Eqs. (68) and (69) can be used to good 
advantage, such as determining the effective room - furniture exchange areas 
(implemented in RMTERM) or to suggest a Gauss-Seidel iterative scheme of 
solving the radiative equations (i.e., a diagonally dominant matrix equation) 


for w andh (implemented in RADFLX). 


4.6.2 Volume Emissivities Due to Soot, H20, and C02 
(Implemented in PLFIRE, WLFIRE, GASUP, and HTAREA) 


The effective absorption coefficient, ger is the absorption coefficient 
in the optically thin limit, or when V is small. To account for a finite 
volume size we extend the mean beam length approximation of Orloff and deRis 


iter to, 


4a V = AL € (70) 
T 
1 oye ey, 
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and L = 3.6 Bey (71) 


where Me is the gas volume, Be is the bounding gas area, and ony is obtained 


and derived from References 13, 18, 19, and 20 to get: 


ep = €. + a - ee é (72) 
Se l - exp (-k_L) ; (73) 
€ = € + € - € : (74) 
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The partial gas pressures within the moving fire plume is approximated as: 
P =MR T,/V (82) 


where He is the volumetric flow rate of the plume via McCaffrey formula [21], 
which can also be found in Reference 2. The plume gases has contributions 
from both the combustion and from the entrained air. The soot extinction 


coefficient, k of the flame is simply Eq. (21). The partial gas pressures 


sf’ 
for the upper and lower gas layer is obtained using the simple gas law. Since 
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the soot cloud exists in the upper gas layer, the method for determining the 
soot extinction coefficient in the upper gas layer is the following. By 
definition the laser soot extinction coefficient is 
-e 
k =C oa (83) 

up ups 
where Eon is the soot mass density obtained from the output of the FAST code 
and mean specific soot extinction area is mass averaged as a function of the 
soot accumulated in the upper gas layer. For radiative heat transfer analysis 
the laser soot extinction coefficient is converted to the mean soot extinction 


coefficient by 


k ~ (Fa), 


i : (84) 
cee n()) 2) 9 tak 


where the averaged wavelength \ is inversely proportional to the upper gas 


layer temperature. 


4.6.3 Evaluation of Surface-gas Exchange Area, oy 


The surface-gas exchange area, TES is the most important of the 
exchange areas in the radiation equations because it dominates the net heat 
transfer rates in Eqs. (63) and (64). That is, the largest temperature 
difference occurs between a flame and a surface, leading to high net heat 
transfer rate. The temperature difference between a surface and another 
surface or between a flame volume to another flame volume is much smaller, and 
implies a much less net heat transfer rate. From Hottel and Sarofim [22] the 


surface-gas exchange area is defined as, 


r(S._) 


x aE fe a 
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and Fic is the distance from a point on A; to a point inside ai Letting ih 


= eo Coe and integrate over the distance “s of the flame volume we obtain, 


T(S. =2> vers.) cos6, cos® 
PES Sg bo ak a cee Ee See 
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As a first approximation we assume mean values for ss Bet: and ae can be 


used when integrating over the flame area ne or we get 


si8, 7 Ja, (Si - 2) €(L Fai dA (88) 
where 
cos®, cos® 
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is the configuration factor between the flame surface, 7 and the surface area 
dA.. Using the reciprocal relationships, F,. dA. = dF -A andAF , = 

i di-y i y-di y ere eS 
A,F, A in Eq. (88) and assuming the surface A; is small enough for r and «€ to 


be constants, we have an approximation to Eq. (88) as, 


es = TOS Ae) ok a) ERNE (90) 


where PS is the configuration factor between the flame surface A and the 
surface area, A,. To make further progress, we address a special case where 
A; = dA. and Eee) = 1 in Eq. (87). Then the surface-gas exchange area, 
Syn! divided by A, is the same as the flame emissivity to an elemental 
surface as defined by Danyan and Tien [23], Modak [24], and Orloff [25]. 
Indeed, these references have numerically computed the flame emissivity from 
an arbitrary flame volume by evaluating the special case of Eq. (87). Orloff 
[25] has thus derived an empirical form for the flame emissivity for normal 


surfaces far from a flame source as: 
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where Er is obtained from Eqs. (72) to (81) and (92), A is the 
ary ary 


projected area of the flame volume, Ne seen by the surface element, and agg 
is the distance from surface element A; to the flame volume centroid. 

Equation (91), however, can be generalized to include the surface orientation 
of A; with respect to a and to include configuration factors and Fick ee ng) as 


Eqs. (88) to (90) seem to suggest. The approximation to Eq. (90) is then, 


(s.g_) 
Sates Ser), Delran bw Sct (93) 
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ie 
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F (94) 
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8, 7 is the angle between the surface element normal vector and the distance 


vector from the centroid of the Me to the surface element. The new projected 
area is thus not calculated directly from the geometry of the flame volume, 


but indirectly from Eq. (94). That is, the configuration factor F, : will be 


directly calculated analytically for a specified flame geometry or 
computationally for a generalized flame shape. Equation (94) is then used to 


* 
calculate an effective projected area, A , which is used in place of A 
i,7 Ly 


in Eq. (92) to obtain e Finally, Eq. (93) is used in place of Eq. (91) 


Thy 


with the objective of obtaining the exchange area eTey for any distance from 


the flame volume or for any optical thickness. 


It is apparent in the optically thick limit, where Er m= rLeU bd 93) 
AS 
reduces to the correct expression for any distance between the flame volume 
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and the surface element. It is also apparent that Eq. (93) will reduce to 
Orloff’s expression for moderate optical thickness when the surface element is 
* 
far from the flame volume, because then we have A = My . For a surface 
i,y WAS 
element located on the flame bounding surface, Eq. (71) is equal to Eq. (92), 
which results in Ps being equal to (0.95 ye Note that if the flame 


ee 
* 
volume is a sphere, then from Eq. (94) we have A on the spherical surface 


bla 
equal to Ay /4.0 for 86. = 0 and F, = 1.0. Thus at least for a spherical 
Y 1,7 1,7 
flame volume, the effective projected area is only 5% in error with that 
expected from the mean beam length theory for the radiative flux from a 
radiating volume of moderate optical thickness to a bounding surface element. 
Thus, for intermediate distances or for non-spherical flame shapes, we can 


expect Eqs. (93) and (94) to also provide good approximations. 


To achieve greater accuracy at the vicinity of the flame base region, we 


have modified 4. = and O. r to correlate with advanced thermal radiation 


calculations for regions around a conical frustrum [2,23-25]. The result is, 


Le mie) cto : (95) 
1,7 1,7 7 
nr. = RUBE Ge ton MO uek + 0.0613 exp¢L/7))| : (96) 
Pe -1 2 
minimum O; a cos (cos (a/2)) : (97) 


a is the distance from the surface element to a point on the flame 
bounding surface, which is on the line connecting the volume centroid and the 
surface element. 1B is described as the distance between the flame base 
center and a "center of radiation" of the flame. r is the flame base radius 


and a is the angle of the conical frustrum surface from the flame base. 


4.6.4 Gas-Gas and Surface-Surface Hottel Exchange Areas 


The evaluation of surface to surface and the volume to volume exchange 
areas are straight forward when making the type of approximations as was done 
for the surface to gas exchange areas. The surface to surface exchange area 


is defined as, 
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assuming a mean value of transmissivity, we obtain 


s,s. = r(2 JA 


isj a tae shy 


If there are several isothermal radiating gas volumes between two 


particular body surfaces Eq. (99) can be generalized to 


=, ad ad ie Or 100 
i) i ? ae | x y aks 
where y. F Ag - F. i (101) 


An alternative to the mean transmissivity, 1, is to use the mean beam 


length approximation, or 


ee = ES, (102) 


In the code the mean transmissivity is converted to ¢e,, via Eq. (86) and 


ik 


Er is then obtained from Eq. (72). Hottel and Sarofim [22] have derived 


similar relationships for s,s,. They also concluded the mean beam length 


concept is useful for ee radiative flux in systems of low optical 
thickness, but not in all systems of high optical thickness. However, when 
the gas in an enclosure is the energy source so that the radiative flux is 
dominated by the gas-surface exchange at high optical thickness, the error in 


s.S, becomes unimportant. 


ij 
The volume to volume exchange area is defined as 
r(S ) 
- k dV k dV : 103 
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Letting Me _ SO, any ae and substituting into Eq. (103) we obtain 
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Let S . be much greater than os or oy and define a mean value for the 


transmissivity and emissivities so that Eq. (104) after integrations and 


substitutions is 


BB metus 
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We see this equation is analogous to Eqs. (90) and (99). Invoking the mean 


beam length approximations, we have the result, 


r(S - (L, + L)) e(L) €(L,) = r(S_ 3 


at te: 4) (2) €(4,) (106) 
Edwards and Balakrishnan [26] have in fact used this approach to obtain 
an approximation for b8,, between remote spherical volumes and again for 
adjacent cubical volumes. Remoteness of the spheres was shown by comparisons 
with the more exact calculation of Hottel and Sarofin [22] for BB. to be only 
as far away as three mean beam lengths, centroid to centroid. Edwards and 
Balakrishnan assumed the remote spheres approximates any two arbitrary remote 
volumes. We generalize this by not limiting the analysis to remote spheres, 
but instead calculating the configuration factor, the gas surface areas, and 
the emissivities between actual volumes. As a consequence, the remoteness 
criteria may be relaxed, and thus a separate By approximation for adjacent 


volumes would be unnecessary. 


By simplifying all the Hottel'’s exchange areas to involve only 
configuration factor, surface areas, and emissivities of both the flame and 
body surfaces has made possible a practical thermal radiation model of finite, 
3-D, gray flames and bodies. The computation effort has been drastically 
reduced, and yet the model is realistic enough to take into account all 


radiation features of a furniture fire. The radiation zonal model has been 
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implemented in the subroutine FLAMES containing the specialized subroutines 
WALLCL, GASUP, HTAREA, RADFLX, SHFAC, RMTERM, and DSHFAC. 


Thus the input to the thermal radiation analysis are the flame and body 
structures and properties and the output are the net thermal radiative powers 
of all gas layers, flames, and surface elements. These output values are fed 
back into the formulae described in this section for a another iteration step 
leading to convergence. The output values are also used by the FAST model for 


the gas dynamics in the roon. 
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5. THE PROGRAM STRUCTURE 


The objectives for developing the program structure of the FAST/FFM were 
the following: (a) implement the equations in Sections 3 and 4 and other 
relevant equations with efficient numerical schemes in the code, (b) obtain a 
very effective coordination of time integrations between FAST and FFM as 
explained in Section 4, (c) develop very few, but effective interface 
subroutines between FAST and FFM, and (d) organize effectively the input 
database of FFM. The accomplishment of these objectives should also make FFM 
more accessible to other versions of room fire models and to any upgrades of 


the FAST model. 


In any complex code involving time integration, such as FAST/FFM, 
selection of an appropriate numerical integration scheme to ensure stability, 
accuracy, and efficiency is very important. Inadequate schemes of numerical 
time integration can lead to subtle but significant truncation errors not 
detectable by the user. To get reasonable predictions, the code may be forced 
to an inordinately small time step that uses too much computer time. Section 
4 also discussed other aspects of numerical time integrations that could lead 
to inaccuracies. The solution of these several problems has resulted in a 
complete revision of the early versions of HEMFAST. The corresponding 
modifications of FFM algorithms were reported in the Interim Report [4] and 
tested successfully. A close comparison between the time integration schemes 


of FAST and FFM guided the development of interfacing subroutines. 


We note that the influence of the radiative heat flux, from the upper 
gas layer (before flashover) on the furniture burnrate, is minor compared to 
the radiative and convective heat flux from the flames to the burning 
elements. Thus most time integration variables in FAST can be allowed to lag 
by a time step with respect to the time integration variables in FFM without 
loss in accuracy or solution stability in FFM predictions. In addition, the 
time integration variables in FAST are quite responsive to the furniture 
burnrate history. These considerations suggest that the FAST modules for time 
integration should be executed directly after calculating the furniture fire 
growth over a time step. A placement of the FAST modules within the Local 
Radiative Heat Flux iteration loop (see Figure 4-1) was avoided because of the 


minor effects of FAST predictions on the furniture burnrate and the species 
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production and because of a desired reduction of computer solution complexity 


and time. 


As another comment, each time-integration variable in FAST is designed 
to be stable and accurate. Likewise, the integration of the virgin surface 
element temperatures in FFM is designed to be well behaved. The formulations 
of the flame spread rate and of the local heat release rate were designed to 
be well behaved and are slowly changing functions of time in most situations, 
so that Euler time integration to estimate the flame front positions and the 
burn history, respectively, should be reasonable. Ensuring stability and 
accuracy in the whole solution is most vexing and has proven to be a roadblock 
to many modelers of room fires or of any other complex physical phenomena. 

The careful attention given to constructing the FAST/FFM functional flowchart 


as shown in Figure 4-1 was done to minimize problems of this nature. 


The program tree shown in Figure 5-1 give the structure for the 
interface between FAST and FFM as well as for the structure of the FFM code 
itself. The dashed boxes specifically enclose the FAST subroutine names, 
while the remaining character strings are the names of the FFM subroutines and 
function routines. Table 5-1 provides a brief description of all the FFM 
subroutines and functions listed in the program tree. The names are arranged 
in the alphabetical order and references are made to their respective logic 
flowchart figures. The logic flowchart figures were arranged into three main 
groupings to follow that of the program tree. The first grouping are the 
flowcharts of the four interface subroutines between FAST and FFM. The next 
grouping of flowchart are arranged around the first phase of the FFM modeling, 
which corresponds to Section 3, Scaling Procedures for Bench Scale Data. This 
flowchart grouping also includes correspondence to subsection 4.1, Furniture 
Geometry Construction. The third grouping of flowcharts correspond to the 
implementation of the rest of the equations in Section 4, The FAST/FFM 


Calculation Systen. 
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TABLE 5-1 
FFM SUBROUTINES AND FUNCTIONS 


SUBROUTINE FLOWCHART 


NAME FIGURE NO. CALLED BY FUNCTION 
BENCHS 5-8 FFMDAT Reads in and processes the bench scale 


BNGEOM 5-21 BNHIST 

BNHIST 5-11 FFMGO 

BNINIT 5-22 BNHIST 

BNSTEP 5-23 BNHIST 

COUNTR -- FFMOUT,FIRPLY, 
FIRPNT , PMERGE, 
SPRDLN , BNGEOM , 
BNINIT, PLFIRE 

FCORNR -- MOCKUP 

FFMDAT 5-2 *NPUTP 

FFMGO 5-3 *SOLVE 

FFMINT 5-4 DSOURCE 

FFMOUT oa) *RESULT 


y Subroutines from FAST code 


data. 


Calculates geometry parameters of any 
existing burning polygonal areas. 


Determines burn histories and 
geometries of burning polygonal areas. 


Initializes burning parameters for 
virgin elements inside burning areas. 


Increments burn history of all burning 
elements over the time step. 


Revolves the indexes by incrementing 
or decrementing the appropriate 
values. 


Locates the corners defining the 
perimeter of a quad. 


Create database for the FAST/FFM 
calculation system. 


Converts FAST input variables to FFM 
variables, executes FF growth modules, 
and outputs data for input to FFMINT. 


Interpolates FFM output variables to 
current values corresponding to FAST 
input variables for SOURCE. 


Furniture fire model output. 
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SUBROUTINE FLOWCHART 
NAME FIGURE NO, CALLED BY FUNCTION 

FIRPNT 5-15 TIFSST Converts ignition point to small 
burning polygon, finds midpoint, 
defines 8 initial flame spread 
vectors. 

FIRPLY 5-16 TIFSST Calculates flame spreading on a 
polygon in a panel. 

FLAMES ere hd FFMGO Generates structures for: flame on 
furniture, upper gas layer, and for 
the room. Calculates combustion 
release rates. 

GASUP 5-29 FLAMES Constructs surface elements on the 
upper gas layer. 

HTAREA 5-30 FLAMES Calculate Hottel’s exchange areas and 
the radiative heat fluxes and powers. 

INTIAL 5-13 TIFSST Starts object burning at given point. 
Also calculates new surface 
temperatures. 

MCORNR -- MOCKUP Transforms panel corners from body 
coordinates to world coordinates. 

MIRROR -- QUADLL Generates a mirror image of triangle. 
Used for left hand triangle. 

MOCKUP 5-6 FFMDAT Generates the chair geometry. Finds 
center point, corners, vectors, and 
true areas. 

MOVE -- MCORNR, Combines action of rotation and 

MOCKUP translation on a position vector in 
grid coordinates. 

PLFIRE 5-26 TOFIRE Generates a merged pool/wall flame 


TABLE 5-1 (continued) 
FFM SUBROUTINES AND FUNCTIONS 


above the given area "I" located in 
the panel "J". 
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TABLE 5-1 (continued) 
FFM SUBROUTINES AND FUNCTIONS 


SUBROUTINE FLOWCHART 


NAME FIGURE NO, CALLED BY FUNCTION 
PMERGE 5-20 TIFSST Merge overlapping flame spread 


polygons, create new spread point when 
distance between points is too large. 


PROCESS 5-9 BENCHS Scales the cone calorimeter data for 
heat, mass, and combustion release 
rates. 

QUADLL 5-7 MOCKUP Describes a quad, returns size of area 
elements, center points, and normal 
vectors. 

RADFLX 5-32 FLAMES Gauss-Seidel iteration for solution of 


Hottel’s zonal equations. 


RECTNG -- QUADLL Returns size, location, and normal 
vector in Body coordinates of area 
elements of a rectangular region. 
Elements are of equal dimension. 


RMTERM eel GASUP Adjusts the Hottel’s exchange areas to 
ensure energy conservation at 
equilibrium temperatures. 


ROTATE -- MOCKUP, Performs rotation of coordinate 
MOVE vectors. 

SHFAC -- GASUP Shape factor subroutine. 

SPRDLN 5-17 FIRPNT Computes flame spread path in given 


direction and the time increments of 
spreading flames. 


SPREAD 5-18 SPRDLN Computes flame spread rate due to 
preignitive heat fluxes for the 
advancing flame front. 


STROUT -- MOCKUP Converts a shape from grid coordinate 
notation to sequential notation. 


ay 


SUBROUTINE 
NAME 


SURTMP 


THMOUT 


TIFSST 


TOFIRE 


TRANLT 


TRIANG 


WALLCL 


WLFIRE 


DSHPAC 


EGAS 


FLOWCHART 
FIGURE NO. 


5-14 


5-24 


5-25 


5-28 


5-127, 


TABLE 5-1 (continued) 
FFM SUBROUTINES AND FUNCTIONS 


CALLED BY 


INTIAL 


FLAMES 


FFMGO 


FLAMES 


MOVE 


QUADLL 


FLAMES 


TOFIRE 


SHFAC 


FFMGO ,HTAREA, 


PLFIRE, PROCESS, 


WLFIRE 


FUNCTION 


Incorporates damping factor to 
temperature predictor so the predicted 
temperature is pointed toward 
equilibriun. 


Determines heat and mass output of a 
single furniture element as a function 
of time since it was ignited. 


Calculates thermal ignitions, flame 
spreadings, and surface temperatures. 


Generates flames on furniture: 
a) merged pool/wall fires, 

b) purely wall fires on either 
vertical or inclined panels. 


Translates vector from body 
coordinates to reference 
coordinates - by amounts of x,y,z. 
(Referenced to the body origin.) 


Returns the given element’s area, 
center point, and normal vectors. All 
areas made up of either rectangles or 
triangles. 


Constructs surface elements for walls 
and ceilings. 


Generates a wall flame on area "I" in 
panel "J". 


Calculates configuration factors 
between two surface geometries. 


Evaluates gas emissivity due to carbon 
dioxide and water. 
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FUNCTION 
NAME 


FA 


FLMVEL 


HC 


FLOWCHART 


FIGURE NO. 


5-19 


TABLE 5-1 (concluded) 
FFM SUBROUTINES AND FUNCTIONS 


CALLED BY 


FFMGO ,HTAREA, 
FFMINT, PLFIRE, 


WLFIRE 


SPREAD 


SURTMP 


FUNCTION 


Coefficients for converting from soot 
volume fraction and gas temperature to 
soot absorption coefficients. 


Calculates flame spread rates. 


Finds convective heat transfer 
coefficients. 
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5.1 Modules Description for Interfacing FFM and FAST 


The FAST/FFM interface is implemented through four FFM subroutines; 
FFMDAT, FFMGO, FFMINT, and FFMOUT. The function of these subroutines was kept 
simple as much as possible for modularity, especially for accomodating the 
updates of FAST. The subroutines FFMDAT (see logic flowchart Figure 5-2) and 
FFMOUT (see logic flowchart Figure 5-5) provide input and output to FFM 
respectively, and they do not share any common block variables with FAST. The 
purpose of the subroutine FFMGO (see Figure 5-3) is to calculate the furniture 
fire behavior over a time step just before DIFEQ and other subroutines called 
by SOLVE calculates the wall and gas layers behavior over the very same time 
step. Another purpose of FFMGO subroutine is to convert FAST variables to FFM 
variables and to save FFM predicted variables to an interpolation common 
block. The subroutine SOLVE also involves sub-time stepping while the FAST 
subroutine DIFEQ (which calls DFE=DSOURC) has the smallest time stepping to 
solve the stiff ODEs of the FAST model. Since DSOURC in turn needs FFM 
predicted variables at the smallest time steps, the call to the FFM subroutine 
FFMINT (see Figure 5-4) will linearly interpolate the FFM predicted variables 
obtained via the interpolation common block. Because FFMINT interpolates at 
the smallest time steps, the subroutine FFMGO can be placed at any level of 
time stepping in the FAST model. Thus during the onset of flashover, the call 
to FFMGO should be in DIFEQ or DSOURC to accomodate the strong coupling 
between the physical processes of the ceiling gas layer with that of the 
furniture fire. Several FAST’s subroutines for defining a fire source, called 
within the DSOURC subroutine, are bypassed to take full advantage of FFM 


capabilities. 


58 


Input: 
Mockup data 
Bench data 
Cone data 

Ignite data 


CALL MOCKUP 
CALL BENCHS 


Output: 
Processed data 
for furniture 

fire model 


MOCKUP: construct 
mockup and ignite 
a spot on mockup 


BENCHS: process 
data from flame 
spread apparatus 
and cone calorimeter 


Figure 5-2 Logic of reading and processing input data to FFM with 
interface subroutine, FFMDAT. 


s}!) 


START 


Input: 
FAST common block variables 
FFM common block variables 
Time and time step 

Old FFM predicted variables 


Reinitialize variables 
Convert FAST variables to FFM variables 
Calculate upper and lower layer emissivities 
TIFSST: thermal ignition, 
CALL TIFsst |flame spreading, and surface 
CALL BNHIST 


temp. updated to current time. 
CALL FLAMES 


Convert FFM predicted variables 
to variables for interpolation in 
FFMINT for input to FAST. 


Print Results 


RETURN 


BNHIST: obtain burn history 
of each burning element. 


FLAMES: calculate radiative 
heat transfer, burnrates, combus- 
tion rates, and surface convective 


heat rates. 


Figure 5-3 Logic of executing furniture fire growth in FAST with the 
interface subroutine, FFMGO. 
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START 


Input: 
Predicted variables from FFMGO 
Input variables from SOURCE 


Time step for interpolation 


Fire source 
in room 


Yes 


No 
Interpolate heat rates to Calculate upper and 
upper and lower gas layers. lower layer emissivities. 
Interpolate flame height. Calculate radiant heat 
transfer. CALL FIRRAD 


Interpolate combustion product 
mass rate to upper gas layer. 


Calculate plume mass rate. 
CALL FIRPLM 


Figure 5-4 Logic for interpolation of furniture fire predictions for FAST 
heat sources in subroutine, FFMINT. 
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Output of furniture fire growth 
from FFMGO in common blocks 


Yes 
RETURN 


Print polygon #, area, 


Panel corner points 
output for plotting 


Sum polygon areas. 


Print fire spread vectors, 
magnitudes, and status. 


<n Yes 


Next polygon 
in panel? 


Fire spread vectors 
output for plotting 


Yes 


No 


Calculate relative burn area of panel. 


Yes 


Ee IPLOT22 ard Output for plot- 


ting burn areas 
and flame foot 
heights on panels 


eC Print surface elements status 


Logic of printing FFM predicted variables in subroutine FFMOUT. 


Next panel 
on mockup 


Figure 5-5 
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5.2 Modules Description For Scaling Bench Data and Initializing the 

Furniture Fire Scenario 

The program tree shows the interface subroutine FFMDAT calling the 
subroutines MOCKUP and BENCHS. The logic flowchart for MOCKUP is shown in 
Figure 5-6 and is the implementation of Section 4.1. Note that MOCKUP was 
designed to be called by FFMDAT, rather than from FFMGO, because the furniture 
geometry and the ignition location only need to be initialized at the 
beginning of the simulation. The subroutine QUADLL (see Figure 5-7) is called 
in turn by MOCKUP to generate surface elements within a trapezoidal panel. 
Several such panels connected together form a furniture mockup. This method 
of furniture construction seem to provide just enough detail to calculate the 
furniture burnrate history and be adaptable to various flame structures while 
at the same time required minimal input from the user to design a furniture 
mockup and still permit practical 3-D thermal radiation calculations. The 
subroutine BENCHS ( see Figure 5-8) basically reads in the bench scale data 
and then processes the cone calorimeter database with a call to the subroutine 
PROCSS (see Figure 5-9). These two subroutines implement Section 3 with the 
manual calibration of the constants. The automatical calibration of the 


constants await incorporation of a couple new processing subroutines. 
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CALL QUADLL QUADLL: generate square- 
like elements in the 


trapezoidal panel and 
provide element positions 
and areas. 


CALL MOVE MOVE: translate panel 
to a location in the roon. 
CALL ROTATE ROTATE: rotate panel to 
the correct orientation. 
CALL STROUT STROUT: reorder panel 
surface elements in 


sequential order and de- 
fine position and normal 
vectors of surface 
elements. 


Define panel parameters. 
CALL FCORNR FCORNR: define corner 
points of panel sequenced 


counterclockwise around 
panel normal vector. 


CALL MCORNR MCORNR: translate and 
rotate corner points to 


panel’s position and 


Yes orientation in the room. 
Next Define connections between adjoining panels. 
panel? 


Yes 


backside 
also? 


Generate opposing panels, 
surface elements, normal 
vectors, positions, areas, 
and define opposing panel 
connections. 


No 


Figure 5-6 Logic of constructing furniture "panels" with subroutine, MOCKUP. 
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START 
Generate left hand triangle 
of quadrilateral. 
CALL TRIANG TRIANG: construct surface 
elements in the triangle 


area. 


Generate right hand triangle 
of quadrilateral. 


CALL TRIANG 


Generate central rectangle 
of quadrilateral. 
CALL RECTNG RECTNG: construct sur- 
face elements in the 


rectangle area. 


CALL MIRROR MIRROR: convert first 
triangle to its mirror 


image. 


Reorder surface element 


positions, normal 
vectors and areas. 


Figure 5-7 Logic diagram for generating a quadrilateral furniture panel with 
subroutine, QUADLL. 
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Input bench scale data 
from a flame spread 
apparatus and the cone 
calorimeter. 


CALL PROCSS 


PROCSS: process the cone 
calorimeter data. 


Figure 5-8 Flowchart for evaluating the bench scale data with subroutine, 
BENCHS. 
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START 


Loop over cone data for 
rm ort rr: different times and 
i irradiance levels. 


| Calculate radiative and 

1 convective heat flux on 

F horizontal samples in cone 
calorimeter due to cone 
heaters and the flame. 


1 

| Calculate soot mass fraction. 
I 

| Calculate quasi-stochiometric 
I heat release flux. 

I 

I 


Integrate quasi-stochiometric heat 
J release flux or fuel mass 

\ release flux to define actual 

1 burn history. 


i 
etl s Pee Dens asd on - End of loop. 


Apply scaling to the heat and 
mass release rates and to the 


burn history for all data points. 


RETURN 


Figure 5-9 Logic flowchart for scaling the cone calorimeter data with 
subroutine, PROCSS. 


5.3 Modules Descriptions For FAST/FFM Calculations 


The implementation of the remaining part of Section 4, FAST/FFM 
Calculation System, is described by the logic flowchart Figures 5-10 to 5-32. 
The prediction of Thermal Ignition, Flame Spread, and Surface Temperature is 
managed by the subroutine TIFSST (see Figure 5-10) called from FFMGO. TIFSST 
in turn call a sequence of subroutines whose logic flowcharts are described in 
Figures 5-13 to 5-20. The implementation of Section 4.3, the prediction of 
local burn history and related geometry parameters, is managed by the 
subroutine BNHIST (see Figure 5-11) called next by FFMGO. BNHIST in turn 
calls the three subroutines, BNGEOM, BNINIT, and BNSTEP to do the detailed 
calculations explained in the logic flowcharts, 5-21 to 5-23. The subroutine 
FLAMES, called last by FFMGO, manages the iterative loop to solve for the 
surface heat fluxes formulated in Sections 4.4 to 4.6. The local convective 
heat flux and the fuel pyrolysis rate for each burning surface element is 
calculated in the subroutine THMOUT (see Figure 5-24). The calculation for 
the chemical heat and mass release rates of the adjoining flames is managed in 
the subroutine TOFIRE (see Figure 5-25). In fact, several merged pool/wall 
flames and wall flames are constructed in this subroutine. The subroutines 
PLFIRE and WLFIRE (see Figures 5-26 and 5-27) do the actual constructing of 
the respective two flames as they are called from TOFIRE. The thermal 
radiation analysis is implemented in the remaining subroutines described in 
Figures 5-28 to 5-32. The remaining surface elements needed for the thermal 
radiation analysis are for the walls and ceiling (WALLCL) and for the gaseous 
layers (GASUP). The Hottel’s exchange areas are calculated in HTAREA and 
subsequently smoothed in RMTERM to preserve energy conservation. RADFLX 


solves the Hottel’s zonal radiation equations with the Gauss-Seidel method. 
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Input: 
Surface element status 


CALL INTIAL INTIAL: obtain any new 
ignition source location 


and the updated element 


temperatures. 
ignition to burning poly- 
gon area. 


CALL FIRPNT 
CALL FIRPLY 


Burning 
polygon? 


burn areas? 
FIRPLY: flame spread 
outward from burning 
surface polygon. 


Select panel 


CALL PMERGE PMERGE: merge any over- 
lapping polygons and re- 
tain nominal distance 

between polygon vertices. 


Yes More No 
panels? 


Figure 5-10 Logic for mockup heating, pyrolysis, ignition, and flame spread 
with subroutine, TIFSST. 
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Vani 


CALL BNGEOM BNGEOM: calculate geometry 
parameters of any existing 


burning polygonal areas. 


CALL BNINIT BNINIT: initialize burning 
parameters for virgin elements 


inside burning areas. 


CALL BNSTEP BNSTEP: integrate scaled burn 
time for each pyrolzing element. 


Figure 5-11 Logic for the burn history of each element with subroutine, 
BNHIST. 
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START 


Input: 
Input data from FAST via FFMGO 
Input data from FFMDAT 


CALL THMOUT THMOUT: get actual heat and mass 
release rates from scaled heat 


and mass release rates using 
surface radiative fluxes 


CALL TOFIRE TOFIRE: construct various pool 
and wall fires over burn areas. 


Define flame temp. and emissivity. 


CALL WALLCL WALLCL: construct wall and ceil- 
ing surface elements. Define 


wall temperature and emissivity. 


CALL GASUP GASUP: construct upper gas layer’s 
surface elements and obtain volume 


and enclosed area of lower gas 
layer. Define temperature and 
emissivity. 


CALL HTAREA HTAREA: calculate Hottel’s ex- 
change areas and ensure energy 


conservation of equilibrium 


Adjust radiative heat fluxes 
and soot volume fractions. 


RADFLX: solve for surface and 
volume radiative power for mock- 
up, room, flames, and gas layers. 


Radiative heat 
rates converge 


Figure 5-12 Logic for calculating the various heat and mass release rates 
of multiple fires with subroutine, FLAMES. 


re 


START 


Input: 
Data from FFMDAT through 
FFM common blocks 


Yes 
User's selection of 
ignition on mockup 


User's ignitiopd 


Loop over 
all elements 


No 


Initialize burning and 
flame spreading parameters 


Non- burning 
element? 


Yes 
CALL SURTMP SURTMP: calculate surface 
temp. for a time increment. 


Yes 


Set temp. to burning value 
and calculate "burning" time. 


above ignition? 


Initialize burning and flame 
spreading parameters. 


Calculate "burning" time. 


Initialize burning parameters. 


Element 
smoldering? 


End of loop 
(ras) 


Figure 5-13 Logic for calculating dynamic thermal status of non-burning 
surface element with subroutine, INTIAL. 


Vz 


START 


INPUT: 
Material data from FFMDAT 
Time step and element } 

radiative heat fluxes 


Calculate convective 
heat transfer coefficient 


Thermally YES Predict surface temperature 
thin OT a with thermal thickness term 
material 


2 


NO 


Calculate convective heat flux 


Predict surface temperature 


with thermal inertia term 


Calculate convective heat flux 


Figure 5-14 Flowchart for surface temperature predictions for either 
thermally thin or thick materials with subroutine, SURTMP. 


START 


INPUT: 
Panel and polygon }# 
Mockup geometry 

Burning parameters 


Initialize flame spread parameters 
Define initial flame spread direction and time step 
‘—--- Loop over flame spread directions 


Call SPRDLN SPRDLN: Calculate a new polygon 
vertex from the increment in the 
flame spread vector. 


Bs Se Se End of loop 


Figure 5-15 Logic for converting an ignition point source to a burning 
polygon area with subroutine, FIRPNT. 
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START 


INPUT: 
Panel and polygon # 
Mockup geometry 

Burning and flame spread parameters 


fi 


Define time step 


Loop over flame spread directions 


Calculate flame spread direction unit vectors 


- - > 


Call SPRDLN SPRDLN: Calculate a new polygon 
vertex from the increment in the 


flame spread vector. 


Rte Se ee End of loop 


Figure 5-16 Logic for calculating flame spreading from a burning polygon in 
a panel with subroutine, FIRPLY. 


yee 


START 


INPUT: 
Flame spread parameters from 
subroutines FIRPNT or FIRPLY 


| 


Adjust flame direction unit vectors if wall fire 


Find intersection with a panel edge 
for a given flame spread direction 


SPREAD: calculate the flame 


Call SPREAD spread rate for a given 
direction 
Compute flame travel distance 


NO Update polygon 
lame reached vertex 
panel edge 


Initialize direction 
unit vectors 


Initialize flame 
spread vectors 


Set polygon vertex to 
panel edge intersection 
and reset other parameters 


Calculate remaining 
time interval into 
adjacent panel 


Adjacent panel 
edge burning 


Figure 5-17 Logic for calculating a new polygon vertex from the increment 
in the flame spread vector with subroutine, SPRDLN. 
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START 


Input: 
Properties for flame spread 

obtained in FFMDAT, parameters 
from SPRDLN. 


ie Loop over elements in the panel. 


Save burning element closest 
| to the flame front position. 


Save virgin element closest 
and in front of the flame 
front position. 


| Calculate averaged virgin 
temperatures and radiative 
| heat flux. 


ar ee et) End of loop 


Calculate thermal radiative heat 
flux at flame front, preignitive 
radiative heat flux, preignitive 


surface temperature, preignitive 
surface length, and effective 
heat transfer coefficient with 
surface orientation. 


FLMSPD=FLMVEL FLMVEL: formulae of 
flame spread rates. 


(une) 
Figure 5-18 Logic for computing flame spread rate distance with subroutine, 
SPREAD. 


START 


flame spread parameters 


Yes 


No 
Thermally 
thin? 


Radiative 
exponential 
distribution? 


Radiative 
exponential 
distribution? 


No 


Standard thermally 
thin formula | 


Standard thermally 
thick formula ~— 


Most general formula 
which also includes 
thermal inertia 


More general formula 
for the radiative 
profile 


(mmm) 


Figure 5-19 Logic of selecting formulae for the flame spread rate in the 
function routine, FLMVEL. 
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START 


Input: flame spread predictions 
from FIRPT and FIRPLY 


Find all flame spread lines 
and identify with nearby 
polygon or flame spread line. 


Find any vertex of a polygon 
inside of another polygon and 
identify these polygons. 


Any 


polygon 
mergings? 


Yes Organize polygons 
to the merged 


polygons. 


No Next merged polygon 


Construct merged 
polygon from original 
polygons by selection 
of outermost vertices. 


Add a polygon vertex between 


adjacent vertices if their 
spacing is too large. 


Remove unnecessary polygon 
vertices on panel edge. 
(nm) 


Figure 5-20 Logic to merge overlapping flame spread polygons and maintain 
reasonable distances between polygon vertices with subroutine, 


PMERGE. 


ve 


START 
INPUT: Flame spread results 


smn) dae ye Loop over all panels 
A ies Loop over burning polygons 
' 

1 

1 Calculate polygon area 

| Calculate polygon perimeter 

' 

; 

Calculate polygon centroid 

1 

\ 

i 

. Calculate polygon hydraulic diameter 

1 

1 

i 

' 

1 

i 

i 

' 

\ 


Determine bottom and top vertex! 
of each burning polygon 


~ 7 = S<e Endiofalcop 


Any polygon 
> 95% of 
panel area 


Reset polygon area to panel 
area and reset other polygon 
derived parameters 


NO 


t------- End of loop 


Figure 5-21 Convert a non-burning element as located inside a burning 
polygon to a burning element with subroutine, BNINIT. 
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START 


Input: flame spread polygons 
and surface elements status 


a Loop over all panels. 


Any 
burning 
polygons? 


Yes 
| 


| Loop over all surface elements in panel. 


Yes Element 
| F, eghci Saparaare burning? 


| No 
| | Loop over all burning polygons in panel. 
| | | Sum up angles between vectors 


connecting non-burning element 
to polygon vertices. 


Sum of 
angles > « 


Yes 


| [est etesene seecus co burning. 
Ce <del 


Figure 5-22 Logic to determine the heat and mass release rate of the 
mockup element as a function of time since its ignition with 
subroutine, THMOUT. 
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START 
INPUT: 


FFM time step, heat and mass 
flux and FFM processed data 


SsoOsn amo eee a4 Loop over all panels 


(eee Loop 


ove 
1 NO 
;— -er eee ce eee we & 


YES 


Setup burning parameters 


Update the scaled times of data sets and 
select a data set with the closest net 
surface heat flux 


r all elements of panel 


Increment the scaled burn time of the 
element by using the corresponding “real® 
time increments to reach FFM time step 


A 
! 3 
o : 


Adjust panel and polygon positions 


lagf- = - = = = = = = = = = = =| End of loop 


Figure 5-23 Logic for creating several pool/wall flames with subroutine, 
TOFIRE. 
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Input: processed data from 
FFMDAT scaled burntime. 


Loop over all burning elements. 


Interpolation for heat of pyrolysis 


Formulae for fuel mass flux, 
convective heat flux, and 
wall flame thickness over 
the burning surface element 


Upper 
horizontal 
urface? 


Formula for fuel mass flux 
and convective heat flux 
over horizontal element 


Interpolation of scaled heat release flux 
and evaluation of net surface heat flux 
Calculate chemical heat release 
flux of burning element. 

End of loop 


Figure 5-24 Logic to determine the heat and mass release rate of the 
mockup element as a function of time since its ignition with 


subroutine, THMOUT. 


83 


Input: mockup geometry 


oe Loop over all panels. 


No Upper 
fMl|- g- - eg ger ekr se Ks horizontal 
surface? 


Yes 
arene Loop over flame spread polygons in panel. 


already with 
a flame? 


No 


CALL PLFIRE PLFIRE: construct a "pool" 
; fire flame over the flame- 


spread polygon and the 


adjacent polygons. 
STs Vance Ew WRT End of loop 


ie 
Si Bae ees Loop over all flame spread polygons. 
1 

1 

\ 

1 Yes 

ne i ee already with 

' a flame? 

i 

i 

; No 

! 

L 


~ | End of loop CALL WLFIRE WLFIRE: construct a wall 
fire flame over the flame 


spread polygon. 


Figure 5-25 Logic for creating several pool/wall flames with subroutine, 
TOFIRE. 
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START 


Obtain surface convective heat rate, 
heat and mass release rate, and the 


mass rate units of ani k. ’ co,, co 
max 


HCL a Ho. H,0 by summing over burning 


elements within horizontal polygon. 


Loop over vertices of horizontal polygon, 
if not already done so. 


Find any adjacent burning vertical 
polygon not already part of a flame 
and not a thermal ignition zone. 


Obtain surface convective heat rate, 
heat and mass release rates, and the 


COCO 
s 2 
max 
HCl, HC, and H,0 for the vertical polygon. 


Sum the rates for the polygons to 
obtain the rates for the attached flame. 


End of loop 


e 
mass rate units of 05: k 


Calculate pool/wall flame height. 


Calculate flame foot lengths for 


vertical burning polygons. 
“(4) 


Heat rate from 
horizontal polygon 
too small? 


Remove vertical Yes 
polygons from 
flame. 


Figure 5-26a First part of logic for constructing a merged pool/wall flame 
over horizontal polygon with subroutine, PLFIRE. 
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ae a a Loop over vertices of horizontal polygon. 


Construct flame base, top, and side elements 


on the strip confined by two adjacent 
vertices, polygon centroid, and flame height. 


Flame area from summing element areas 
fk ee ee ee 
End of loop 
Calculate flame volume. 


Formulae for: flame mean beam length and 
emissivity, plume volumetric rate at solid 
flame height, soot mass fraction and mass 


flow rate, partial pressures of co. and H,0, 


and flame total emissivity coefficient. 


Figure 5-26b Second part of logic for constructing a merged pool/wall 
flame over horizontal polygon with subroutine, PLFIRE. 
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START 


Obtain surface convective heat rate, heat 
and mass release rate, and the mass rate 


units of o. k, , €O,, CO, HCl, HC, and H,0 


2 
max 


by summing over burning elements within 
inclined polygon. 


Calculate "solid" height of wall flame. 


-— -— — — | Loop over burning elements. 


Construct surface elements at flame base 
and at flame thickness. 
Sum for the flame area and volume. 
segs 2 oe 
End of loop 


Obtain surface elements for flame foot 
extending to flame height, also the sums 
for the flame area and volume. 


Formulae for: flame mean beam length and 
emissivity, plume volumetric rate at solid 
flame height, soot mass fraction and mass flow 
rate, partial pressures of co, and H,0, and 


flame total emissivity coefficient. 


Figure 5-27 Logic for generating a wall flame on a inclined burning 
polygon with subroutine, WLFIRE. 
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Data input: room 
geometry from FAST. 


Surface elements as defined by temperature, 
emissivities, positions, areas, and normal 
vectors are for: one ceiling element, four 

upper wall elements, four lower wall elements, 

and one floor element. 


Figure 5-28 Logic to create surface elements on the walls and the ceiling 
with subroutine, WALLCL. 
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START 


Input: geometry of gas 
layers from FAST 


Construct elements defined with areas, 
positions, and normal vectors for: 
one ceiling sided element, four wall 
sided elements, five bottom sided 
elements including the square one 
above the fire source. 


Flame 
penetrating 
upper gas layer? 


Yes 


Construct indentation of gas layer 
with surface elements for: dis- 
placed square element above fire, 
and four wall indented elements. 


(sem) 


Figure 5-29 Logic for generating surface elements for the upper gas layer 
with subroutine, GASUP. 
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START 


Set radiative properties of gas zones. 
Ignore upper or lower gas layers they are radiatively clear. 
ah ee or Double loops over material surface elements 


Set up elements for viewfactor with blocking. 


CALL SHFAC 


SHFAC: calculate the view 
factor between selected 
surface elements including 
effects of blocking (pub- 
lished program). 


Calculate surface-surface exchange areas. 


Loop over gas zones (flames and gas layers). 


Include gas zone elements as blocking. 
CALL SHFAC 


eS 
vt aa 


CALL SHFAC 
Calculate volume-volume exchange area. 


Adjust volume-volume exchange area due to lower gas attenuation. 


= Soe eB o eae ee End of loops cs 


Figure 5-30a First part of logic for calculating the Hottel’s exchange 
areas for all regions in the room with subroutine, HTAREA. 


90 


Loop over material surface elements. 


| Set up material viewfactor elenents. | 


Loop over gas zones. 


| Set up gas zone viewfactor elements. | 
| CALL sHFAC | 


Calculate volume-surface exchange areas. 
----co- Loop over gas zones. 
Include gas zone elements as blocking. 
CALL SHFAC 


Adjust volume-surface exchange area 
due to attenuated blockings. 


j= --t—-+------- -- = End of loops 


CALL RMTERM RMTERM: calculate any remain- 
ing Hottel'’s exchange areas 

and smooth out exchange areas’ 
values for energy conservation 


Hottel’s exchange 
Figure 5-30b Second part of logic for calculating the 
oe areas for all regions in the room with subroutine, HTAREA. 


es 


START 
Calculate smoothing factors for all 
gas zones to conserve energy. 


Combine smoothing factors to adjust 
the volume-surface exchange areas and 
the volume-volume exchange areas. 


No Smoothing 
factors 
converged? 


Yes 


Calculate smoothing factors for all 
material surface elements to conserve 
energy. 


Combine smoothing factors to adjust the 
surface-surface exchange areas. 


Smoothing 
factors 
converged? 


Yes 


Figure 5-31 Logic to smooth out the values of Hottel’s exchange areas to 
satisfy energy conservation with subroutine, RMTERM. 
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Solve for surface irradiances of all 
material surface elements by Gauss- 
Siedel matrix method. 


Solve for volume irradiances of all 
active gas zones by Gauss-Siedel 
matrix method. 


Irradiances 
converged? 


Calculate radiances and net radiative 
powers of all material surface elements 
and of all active gas zones. 


Figure 5-32 Logic to solve the Hottel’s zonal radiation equations by the 
Gauss-Siedel method with subroutine, RADFLX. 
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6. DESCRIPTION OF DATA FILES USED BY FAST/FFM 


The management of the FFM database requires stringing together separate 
input data files. One input file provides definition of the furniture shape 
and position; another file provides the complete cone calorimeter database for 
different materials; and yet another provides data input for thermal ignition 
and flame spread. This data management required modification to the FAST 
version 18 database file so that the user can select the different input data 


files for stringing together. 


Dr. Walter Jones of the Center of Fire Research has made the 
accomodating changes to the FAST database file and to the FAST input 
algorithms while converting to FAST version 18 in the current HEMFAST version. 
Accomodating changes in a few FFM subroutines were also made by us to make the 
modified code, now called FAST/FFM, also workable for a wall fire scenario in 


a kitchen. This guide documents these changes to the input data files. 


6.1 FAST Database File 


The reader is advised to be familiar with Section 6, Description of the 
Data File Used by FAST in Reference 1. The basic modification to this 
documentation for the FAST/FFM code are additions to subsection 6.7, Fire 
Specifications of Reference 1. Furniture fire is specified by placing the 
value 5 to the right of LFBT in the FAST input data. Using any other value 
will default the FAST/FFM code effectively back to the FAST version 18. On 
the other hand, if the furniture fire source is specified, the user is 


required to include the following input information: 


Label Parameter Comments Units 


LFBO Room of Room of fire origin is the room number in 
Fire Origin which the fire originates. Default is l. 


LFBT Fire Type Set to 5 for furniture fire 
LFPOS Fire Corresponds to location of center point 
Position of furniture mockup in the room 
1 Center of the room 
2 Corner of the room 
3 Center of the wall, but not at corner 
CHEMI (6) Chemical kinetics and miscellaneous prameters 
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Molar Use default value 


weight 

Relative The initial relative humidity in the system v4 
Humidity 

Limiting The limit on the ratio of oxygen to other % 
oxygen gases in the system below which a flame will 

index not burn. 

Heat of Use default value because information is J/kg 


Combustion ignored in the furniture fire model 


Initial Set to ambient temperature K 
Fuel 

Temperature 

Gaseous Minimum temperature for ignition of the fuel K 
Ignition as it flows from a compartment through a vent 


Temperature into another compartment. 


LFMAX Number of Set to 1 for furniture fire model 
Intervals 
FTIME Time For the furniture fire this is redefined to be s 
Interval the synchronized time step of FFM and FAST. 
(LFMAX) 
FMASS Mass Loss Set to zeroes for the furniture fire case. kg/s 
Rate These values are not used in FFM. 
(LFMAX+1) 
FHIGH Fuel Height Not used for the furniture fire model m 
FQDOT Heat release Not used for the furniture fire model W 
rate 
FAREA Fuel area Not used for the furniture fire model m*m 
MOCKP (N) The label MOCKP indicates that the names of 


mockup geometry files on this line describe the 
furniture fire source. 


BENCH (N) The label BENCH indicates the names of the 
bench scale thermal, ignition, and flame spread 
properties file is on this line. 


CONEF (N) The label CONEF indicates that the names of the 
cone calorimeter database file at a single 
cone irradiance level and for the same material 
as referenced by BENCH and MOCKP is on this line. 
Note CONEF can be repeated to access a cone calor- 
imeter database file for a different cone irradiance. 


fs 


LFPOS 1 

CHEMI 0.0 50.0 10.0 18000000 300 
LFMAX 1 

FMASS 0 0 

FTIME 20.0 

MOCKP MOCKUP1.DAT 

BENCH FRPUHO.DAT 

CONEF CHOLCA. DAT 


In the example, the furniture mockup fire (LFBT 1) originates in room 
number 1 (LFBO 1) in the center of the room (LFPOS). The initial relative 
humidity is 50%, the molecular weight is 16 (zero is not allowed, so the 
default is used) and the limiting oxygen index is 10%. The heat of combustion 
of 18000000 J/kg and the fuel temperature of 300°K are not used in the 
furniture fire model. The synchronized time step is 20 seconds (FTIME 20.0). 
Note that values with LFMAX and with FMASS are required to access the value 
with FTIME. The furniture geometry data is in MOCKUP.DAT file, some bench 
scale data in in FRPUHO.DAT file, and the cone calorimeter data is in the 


CHOLCA.DAT file. 


6.2 MOCKP Database File 


The data required by the furniture fire model to generate the furniture 
mockup geometry and define the ignition location is the user supplied MOCKP 
Database File. The subroutine MOCKUP reads the data in free format and in the 
sequential order for the variables defined and listed in Table 6-1. These 
coded variables are self explanatory and can be found directly in the 
subroutine MOCKUP. Note that only up to four panels are allowed in the 
current version of MOCKUP subroutine. This subroutine is being revised to 
permit more panels and with variation of the furniture material. A more 


convenient method to define an ignition location is also being devised. 
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VARIABLE 


NAME 


IRMNUM (4) 


ELTHCK (4) 


SEATF 


SEATB 


SEATLN 


ALPHAS 


DELTAS 


RIGHTT 


RIGHTB 


RIGHTH 


ALPHAR 


DELTAR 


BACKT 


BACKB 


BACKH 


ALPHAB 


DELTAB 


TLEFT 


BLEFT 


HLEFT 


TABLE 6-1 


INPUT VARIABLES OF MOCKUP DATABASE FILE 


TYPE 


FREE FORMATTED 


Integer array 


Real array 
Real 
Real 


Real 


Real 
Real 
Real 
Real 


Real 


Real 
Real 
Real 
Real 


Real 


Real 
Real 
Real 
Real 


Real 


DESCRIPTION 
Defines the room # for placing the panel. 
A value of zero removes panel from mockup. 
The thickness of the mockup panel (m) 
Long base length of seat trapezoidal panel (m) 
Short base length of seat trapezoidal panel (m) 


Length between the short and long parallel 
bases of seat trapezoidal panel (m) 


Left lower angle of seat trapezoid ( < 90° ) 
Referenced area for seat element areas (m*m) 
Long base length of right trapezoidal panel (m) 
Short base length of right trapezoidal panel(m) 


Length between the short and long parallel 
bases of right trapezoidal panel (m) 


Left lower angle of right trapezoid ( < 90° ) 
Referenced area for right element areas (m*m) 
Long base length of back trapezoidal panel (m) 
Short base length of back trapezoidal panel(m) 


Length between the short and long parallel 
bases of back trapezoidal panel (m) 


Left lower angle of back trapezoid ( < 90° ) 
Referenced area for back element areas (m*m) 
Long base length of left trapezoidal panel (m) 
Short base length of left trapezoidal panel(m) 
Length between the short and long parallel 


bases of left trapezoidal panel (m) 
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TABLE 6-1 (concluded) 


INPUT VARIABLES OF MOCKUP DATABASE FILE 


VARIABLE TYPE 
NAME FREE FORMATTED DESCRIPTION 

ALPHAL Real Left lower angle of left trapezoid ( < 90° ) 

DELTAL Real Referenced area for left element areas (m*m) 

X(4) Real array Length to translate panel in X direction (m) 

Y(4) Real array Length to translate panel in Y direction (m) 

Z(4) Real array Length to translate panel in Z direction (m) 

PHI (4) Real array Angle of rotation around new Z axis measured 
clockwise from new x axis (degrees) 

THETA(4) Real array Angle of rotation around new X axis (having 
been rotated by PHI() ) measured clockwise 
from the Z axis (degrees) 

PSEC4) Real array Angle of rotation around Z axis (body) measured 
clockwise from X axis (degrees) 

IADJ(I,J,K) Integer array Panel # connected to panel I along the segment 
between corner points J and K of panel I 

IOPPOS Integer A value of 1 specifies double sided panels 
A value of 0 specifies single sided panels 

THICKM Real Thickness between double sided panels (m) 

NFIRIG Integer Number of elements in which to ignite 
simultaneously 

IFIRIG(9) Integer The indexes of surface element in which to 


ignite the furniture fire 
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6.3 BENCH Database File 


The data for the material thermal, ignition, and flame spread parameters 
as well as some miscellaneous parameters is the user supplied BENCH Database 
File. The subroutine BENCHS reads the data in free format and in the 
sequential order for the variables defined and listed in Table 6-2. Most of 
the coded variables are self-explanatory and all can be found directly in the 
subroutine BENCHS. Eventually about half the coded variables will be 
hardwired and a few others will be automatically calibrated. In the present 
database, the variables likely to have the most difficulty with the user are 
the following: FOOTLT, QZERO, and SKMAX. Their equivalent mathematical 
and k. 


symbols are 2 and were introduced in Section 3, the scaling 


f? 1p: ,max 

procedures for bench scale data. The reason for the difficulty is that these 
parameters need to be determined in an iterative process and are critical to 
the success of scaling the bench level data. Thus a hierarchy of FFM output 
was developed with the parameter IPLOT (see Table 6-2) to assist the user in 
the manual iterative process of calibrating these three constants. The FFM 


output are fully explained in Section 7, Data Management of Model Output. 


6.4 CONEF Database File 


The cone calorimeter data is supplied by the cone calorimeter apparatus. 
Sometimes, as the cone calorimeter is upgraded for some reason, the format and 
the content of the cone calorimeter database are also changed. To adapt to 
the various formats we have devised a user supplied data to be inserted as the 
first line of the cone calorimeter database. This user supplied data consists 
of the first six variables listed in Table 6-3. At the present we have 
identified three formats for the cone database. To decide on a format for a 
particular cone database the user will have to inspect the database and to 
inspect the read statements in the subroutine BENCHS of the FFM code. Thus 
the user is required to select values for IFRMAT, TSTART, and TEND while the 
values for QR() , IQUAL, and BNTHCK() should be contained in the cone database 
and are required for the furniture fire model. The rest of the variables in 


Table 6-3 are self-explanatory. 
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TABLE 6-2 


INPUT VARIABLES OF BENCH DATABASE FILE 


VARIABLE TYPE 
NAME FREE FORMATTED DESCRIPTION 
EMISIV Real Material surface emisivity 
AIRAK Real Critical heat release rate ratio of short pool 
fire to merged pool/wall fire for flame shapes 
RSTOCH Real Stochiometric ratio of fuel to air 
RMTMPB Real Room temperature at beginning time (°K) 
THRMLT Real Thermal thickness of material (kJ /°Km?) 
THRMIA Real Thermal inertia of material (kW/m? °K)? 5 
PDEL Real Characteristic length of exponentially decaying 
surface radiative flux of adjoining flame (m) 
TIG Real Piloted ignition temperature of material (°K) 
QCFOOT Real Surface convective heat flux of anchored 
diffusion flame at the flame front (kW/m*m) 
FOOTLT Real Effective flame foot length in opposed flow (m) 
QZERO Real Constant surface heat flux parameter for 
scaling heat and mass release rates (kW/m*m) 
SKMAX Real Maximum flaming soot absorption coeffcient(1/m) 
XLAM Real Laser wavelength of the soot optical probe (um) 
SIGE Real Specific soot extinction area (m*m/kg) 
ISTOCH Integer If 0, scale the chemical heat release flux 
If 1, scale the quasi-stochiometric heat flux 
LHUST Integer If 0, the scaled time is the real burn time 
If 1, the scaled time is the heat release/area 
If 2, the scaled time is the mass release/area 
CQ Real First scaling parameter for time stretching 
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TABLE 6-2 (concluded) 


INPUT VARIABLES OF BENCH DATABASE FILE 


VARIABLE TYPE 
NAME FREE FORMATTED DESCRIPTION 

DQ Real Second scaling parameter for time stretching 

DENAIR Real Density of the air (kg/m*m*m) 

HCPAIR Real Height ratios of solid flame to visible flame 

TBURN Real Surface temperature of burning surface (°K) 

GRAV Real Gravitational acceleration (m/s*s) 

SCAT Real Scattering coefficient of smoke 

FLMTMP Real Thermal radiation temperature of flames (°K) 

IRMTRM Real If 0, first mode of thermal radiation solution 
If 1, second mode of thermal radiation solution 

SIGMAG Real Critical value of surface gradient to determine 
pool-type flames versus wall flames. 

PCO2ST Real Stochiometric pressure of CO2 in the flame (not 
used if cone calorimeter data includes C02) 

PHZ0sT1 Real Stochiometric pressure of H20 in the flame (not 
used if cone calorimeter data includes H20) 

IPLOT Integer If 0, only FAST output sent to unit 6 


If 1, also send FFM output to unit 6 

If 2, also send FFM output to units 10,11,12 to 
create files for comparison with data 

If 3, also send FFM output to unit 13 to create 
input file for 3-D flame spread graphics 

If 4, also send FFM output to units 14 and 15 
to assist in calibrating bench scale 
parameters 
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VARIABLE 


NAME 


IFRMAT 


QR() 
TSTART 
TEND 


IQUAL 


BNTHCK() 


QSUM(, ) 
HCOMB (, ) 
W(,) 
FAIR(,) 
FCO2(,) 
FCO(,) 
FH20(, ) 
FHC(, ) 
FHCL(, ) 


U(,) 


TABLE 6-3 


INPUT VARIABLES OF CONE DATA FILE 


TYPE 


FREE FORMATTED 


Integer 


Real array 
Real 
Real 


Integer 


Real array 


Real array 
Real array 
Real array 
Real array 
Real array 
Real array 
Real array 
Real array 
Real array 
Real array 
Real array 


Real array 


DESCRIPTION 
Set to 1,2, or 3 to select format of the cone 
calorimeter data being read in 
Irradiance flux of the cone heater (kW/m*m) 
Time of material ignition in cone calorimeter 
Time of material burnout or extinquishment(s) 


Index describing orientation of cone sample 
(O for horizontal and 1 for vertical) 


Effective material thickness equivalent to the 
simulated panel thickness on mockup (m) 


Time (s) 

Chemical heat release flux (kW/m*m) 

Chemical heat release/area (MJ/m*m) 

Chemical heat of combustion (MJ/Kg) 

Material mass loss rate (g/s) 

Soot extinction area (m2 /kg) 

Carbon dioxide mass fraction of material loss 
Carbon monoxide mass fraction of material loss 
Water mass fraction of material mass loss 
Hydrocarbon mass fraction of material loss 
Hydrocloride mass fraction of material loss 


Total mass flow rate in the duct (Kg/s) 
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7. STRUCTURE OF FFM MODEL OUTPUT 


The parameter IPLOT is the input variable in the BENCH database file in 
which the user specifies the FFM output. If IPLOT is 0, then only the usual 
FAST tabular output is sent to the Unit 6 file. But if IPLOT is 1, then the 
FFM tabular output is intersperse with the FAST tabular output. The FFM 
tabular output typically provide predicted values for the thermal heating and 
flame spreading behavior, for the surface element burning behavior, and for 
the convective and thermal radiative heat rates. An example of the combined 


FAST and FFM tabular output is shown in Appendix B. 


Whenever IPLOT is set equal or greater than 2, then specialized output 
files are created for the sole purpose of plotting the model results. Some 
plotting examples are given in Appendix C. For IPLOT equal to 2, three types 
of validation output files are created and are in turn used by a plotting 
program to compare the model with the full scale data. The first validation 
file (to Unit 10) is simply the tabular output of the FFM predicted variables 
needed by FAST, such as furniture heat and mass release rates and the specie 
production fractions (see Table 7-1). The second validation file (to Unit 11) 
is just the tabular output of FAST predicted variables needed by FFM, such as 
the upper gas layer thickness, the lower and upper gas layer temperatures, 
etc (see Table 7-2). The third validation file (to Unit 12) is the tabular 
output of FFM predictions of the fractional burn areas of the panels (see 
Table 7-3). The comparison of the fractional burn areas with the flame spread 
data from the flame spread apparatus provide a means to calibrate the flame 
foot size, Ae. 

A three-dimensional flame spreading on a furniture is graphically 
presented with a 3-dimensional plotting program when it reads in the graphic 
file (ie. Unit 13) generated from setting IPLOT to 3 (see Table 7-4). The 
validation files are still generated as described in the previous paragraph. 
Finally, to assist in the manual calibration of the parameters qo and eae: 
the plots shown in Figures 3-3 to 3-6 were created with the plotting program 
using the bench scaling files (i.e., Units 14 and 15) generated from setting 
IPLOT to 4 (see Table 7-5). Our plotting programs are user friendly and only 
require simple tabular data input. However, these plotting programs contain 


licensed software which we would not transport. It is our understanding that 
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Dr. Walter Jones has access to equivalent plotting programs in order to plot 


the same figures located in Appendix C. 


104 


TABLE 7-1 


OUTPUT VARIABLES TO UNIT 10 FROM SUBROUTINE FFMGO 


VARIABLE TYPE 
NAME FREE FORMATTED DESCRIPTION 
CTIME Real Current simulation time (s) 
RLHEAT (1) Real Heat release rate of mockup in room I (kW) 
QCFURN (1) Real Convective heat rate on mockup in room I (kW) 
RLMASS (1) Real Fuel release rate of mockup in room I (kg/s) 
F2 Real O, mass fraction of fuel as depletion (%) 
F3 Real 06, mass fraction of fuel as products (2%) 
F4 Real CO mass fraction of fuel as product (%) 
F6 Real HC2£ mass fraction of fuel as product (4) 
F7 Real HC mass fraction of fuel as product (%) 
F8 Real H,0 mass fraction of fuel as product (4%) 
F9 Real C mass fraction of fuel as product (2%) 
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TABLE 7-2 


OUTPUT VARIABLES TO UNIT 11 FROM SUBROUTINE FFMGO 


VARIABLE TYPE 
NAME FREE FORMATTED DESCRIPTION 
CTIME Real Current time (s) 
QGAS (1,1) Real Net radiative power from mockup fire (kW) 
QGAS (2,1) Real Net radiative power from upper gas layer (kW) 
QGAS (3,1) Real Net radiative power from lower gas layer (kW) 
TG(1) Real Upper gas layer temperature (°K) 
TG(2) Real Lower gas layer temperature (°K) 
TW(1) Real Ceiling temperature (°K) 
TW(2) Real Upper wall temperature (°K) 
TW(3) Real Lower wall temperature (°K) 
TW(4) Real Floor temperature (°K) 
SIGAVG*F9/100. Real Soot extinction area Gare) 
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TABLE 7-3 


OUTPUT VARIABLES TO UNIT 12 FROM SUBROUTINE FFMOUT 


VARIABLE TYPE 

NAME FREE FORMATTED DESCRIPTION 
CTIME Real Current time (s) 
AREA(1) Real Fractional burn area of seat panel 
AREA (2) Real Fractional burn area of left arm panel 
AREA (3) Real Fractional burn area of back rest panel 
AREA (4) Real Fractional burn area of right arm panel 
FOOTHT (1,1) Real Wall fire height of seat panel 
FOOTHT (2,1) Real Wall fire height of left arm panel 
FOOTHT (3,1) Real Wall fire height of back seat panel 
FOOTHT (4,1) Real Wall fire height of right arm panel 
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TABLE 7-4 


OUTPUT VARIABLES TO UNIT 13 FROM SUBROUTINE FFMOUT 


VARIABLE 
NAME 

NUMPLT 

NUMPTS 


CORNER (II ,JJ ,KK) 


NPOINT 


SPRDSVC(II,JJ,KK) 


TYPE 

FREE FORMATTED DESCRIPTION 

Integer Number of curves 

Integer Number of points for curve of panel 
Real Corner 3-D points of panel II at 

corner JJ 

Integer Number of points in a burning polygon 
Real Vertices of a burning polygon JJ in 


panel II 
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TABLE 7-5 
OUTPUT VARIABLES TO UNITS 14 AND 15 FROM SUBROUTINE PROCSS 


VARIABLE TYPE 
NAME FREE FORMATTED DESCRIPTION 

NSET Integer Number of separate cone databases at 
different irradiances (Units 14 and 15) 

NPT(I) Integer Number of points in a single cone database 
(Units 14 and 15) 

25SEC); 1) Real Scaled burn time aver of point J in 
database I (Units 14 and 15) 

Qsc(J,1) Real Scaled heat release flux of point J in 
database (Unit 14) 

WSC(J,I1) Real Scaled fuel release flux of point J in 


database I (Unit 15) (g/J) 
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8. NOMENCLATURE 


OTs bounding gas area aa 
A, flame surface area eas) 
A, surface element area ae 
Ap projected area of gaseous volume, Cay 
a 

ae effective absorption coefficient (1/m) 

oe mass concentration of soot in upper gas layer (keyun) 

Cc material heat capacity (J/g°K) 

aoe air heat capacity (J/g°K) 

D hydraulic diameter of the flame base (m) 

d 2nd calibration constant for burn history in Eq. (59) 

e lst calibration constant for burn history in Eq. (59) 

a blackbody surface radiative heat flux (oT 7) (kW/m?) 

e elemental volume blackbody radiative heat flux (or ) (kW/m?) 
Fs Radiation factor define by Bard & Pagni [13] 

Ps carbon monoxide mass fraction of fuel 

F, mass fraction of combustion product, j 

Be eee soot mass fraction of fuel from cone calorimeter 

G time stretching parameter in Eq. (59) 

Gry Grashof number defined by Eq. (46) 

4 gravitational acceleration (u/s7) 

BB, exchange area between the volume Me and the volume Fes ne 
Se) exchange area between the volume Ke and the surface A, no 
H. heat of reaction of carbon (kJ/kg) 

Hoo heat of reaction of carbon monoxide (kJ/kg) 
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rye 


effective heat of pyrolysis (inverse of mn’) (kJ/kg) 


linearized heat transfer coefficient (kW/m2°K) 


convective heat transfer coefficient (kW/m? °K) 

chemical heat release flux from cone calorimeter Gila 
incoming surface radiative heat flux acueuan 

incoming elemental volume radiative heat flux Gana) 


effective heat transfer coefficient in a pool fire (kW/m?°K) 


radiative heat transfer coefficient (kW/m?°K) 

scaled heat release flux (-) 

quasi-stoichoimetric heat release flux (kW/m*) 

material thermal conductivity (kW/m°K) 

air thermal conductivity (kW/m°K) 

soot absorption coefficient in flame Gis) 

maximum soot absorption coefficient for the material ES 
mean soot extinction coefficient in upper gas layer Er 

laser soot extinction coefficient in upper gas layer (aes) 

gas zone extinction coefficient Aes 

mean beam length (m) 

convective contribution to the effective heat of pyrolysis (kJ/kg) 
conductive thermal length (m) 

distance from surface element A; to gaseous volume centroid (m) 
plume radiative thermal length (m) 

fuel mass release rate (kg/s) 

flowing gaseous zone mass rate (kg/s) 

combustion product mass rate (kg/s) 


upward soot mass flow rate (kg/s) 
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mass flow rate of fire plume (kg/s) 

scaled fuel release flux (kg/kJ) 

material mass loss flux from cone calorimeter Ge fats) 
ambient air pressure (atm) 

partial vapor pressure of carbon dixoide (atm) 

partial vapor pressure of water vapor (atm) 

ambient Prandtl’s number (0.72) 

heat release rate of fire (kW) 

total heat release rate of flame (kW) 


heat release rate from horizontal surfaces (kW) 
heat release rate from inclined surfaces (kW) 


surface net radiative source power (kW) 

volume radiative source power. (kW) 

surface convective heat flux (kW/m) 

gaseous surface conductive heat flux at flame front Cyne) 
scaled heat flux constant parameter (kW/m?) 

net surface radiative heat flux Cew/ae) 

surface radiative heat flux emitted from material (kW/m) 
surface radiant fluxes from external sources (ews) 
surface radiative heat flux from attached flame (kW/m?) 
surface radiative heat flux from cone heater (kW/m) 

net surface heat flux (KW/m") 

air gas constant 

a partial gas constant 

hydraulic radius of a flame base (m) 


stoichiometric fuel/air ratio 
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S. 
iy 


distance from a point on A; to a point inside es (m) 
surface distance from the flame front (m) 

preheated distance (m) 

exchange area between surface AS and the surface A, nan 
scaled time (effective heat release per area) eae 
ambient temperature (°K) 

temperature of burning surface (°K) 

preheated surface temperature (°K) 

surface equilibrium temperature (°K) 

temperature in combustion zone of diffusion flame (°K) 
gaseous zone temperature (°K) 

surface ignition temperature (°K) 

material surface temperature (°K) 

initial surface temperature (°K) 

time (seconds) 

opposed air speed at flame front (m/s) 

flame front quasi-steady speed (m/s) 

fire plume volumetric flow rate Ee) 

isothermal finite volume (a) 

outgoing surface radiative heat flux (kW/m?) 

outgoing elemental volume radiative heat flux (kW/m) 
surface distance in flaming boundary layer (m) 

wall pyrolysis zone width (m) 

dimensionless fuel release flux in Eq. (34) 

merged pool/wall fire height (m) 


averaged visible pool fire height (m) 


BY 


solid flame height (m) 


averaged visible wall fire height (m) 

wall pyrolysis zone length (m) 

defined by Eq. (41) 

effective material thermal thickness (m) 

wall flame laminar boundary layer thickenss (m) 
wall flame turbulent boundary layer thickness (m) 
surface emissivity of burning surface 

carbon dioxide gaseous emissivity 

emissivity correction due to C02 and H20 band overlap 
flame emissivity 

gaseous emissivity due to co, and HO 

water vapor gaseous emissivity 

surface element emissivity 

gaseous zone or flaming soot emissivity 

total gaseous zone emissivity 

defined by Eq. (44) 

inclined surface angle of wall fire 

wavelength (m) 

defined by Eq. (42) 

ambient viscosity (kg/ms) 

ambient kinematic viscosity Coley 

volume scattering coefficient 

material density (e/a) 

air density (ke/m>) 


surface element diffused reflectivity 


114 


Stefan-Boltzmann constant (kW/m? °K‘) 

specific soot extinction area Gic /kB) 

soot extinction area per unit soot mass (ms/EE) 

defined by Eq. (43) 

angle between the surface element normal vector and the distance 


vector from the centroid of the ve to the surface element (rad) 


be) 


10. 


ie 


2. 


34 
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APPENDIX A: SAMPLE INPUTS TO FAST/FFM 


VERSN 
TIMES 
TAMB 

HI/F 

WIDTH 
DEPTH 
HEIGH 
HVENT 
CEILI 
WALLS 
FLOOR 
CHEMI 
LFBO 

LFBT 

LFPOS 
FTIME 
THRMF 
BENCH 
MOCKP 
CONEF 


TABLE A-1: WALL.DAT INPUT FILE TO FAST 


18 -1 wall burn test for a single compartment - restaurant 
40 HOR FO0REO 
300. 101300. 0. 
0.00 
20.0 
Poe0 
1.8 
lin 2selge5228 20 128087S.0F07g87 0900 
GYPSUM 
GYPSUM 
CONCRETE 
rGr 0. 10.0 14000000. 300. 500. 
1 
5 
Hi 
108 
THERMAL. TPF 
FRPUHO . DAT 
MOCKUP1 . DAT 
CHOLCA . DAT 
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TABLE A-2: FRPUHO.DAT INPUT FILE TO FFM AS PRESCRIBED IN WALL. DAT 


1.0 0.0 Cyt oe 007 

O7720° 0200128760 .015 71183 2255080133 

30.0 tesn2e 06652895300. Selma ta0 <1, “C120 
1.2845 0.4 700.0 9.82 0. OkseE25C..@ 0 

0.25 OMPGGgeeO, 1215) 4 


19 


MOCKUP1.DAT INPUT FILE TO FFM AS PRECRIBED IN WALL.DAT 


TABLE A-3: 


0.0413444 


0.0413444 
0.0413444 
0.0413444 


Added 


tooo 


ONO SO 


NON O 


OoOmMo + 


ooo oO 


SY 6.0) (2) S&S 


OM © 


On OO 


NOOO 


ot Ooo 


Cab) Sp tS 


OA ON 


MOCO 


oOoO0oC oO 


oHOooo 


on OM 
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TABLE A-4: CHOLCA.DAT INPUT FILE IN FFM AS PRESCRIBED IN WALL. DAT 


20725°- 240% 474R Gah 0.05 
HEAVY OLEFIN/CA HORIZONTAL 25KW/M2 27.2s 11/18/82-244 (CAL.CONST. 


TIME Q-DOT SUM HEAT OF MASS ~~ AIR/ co2 CO 
Q-DOT COMB LOSS SAMPLE  FRACT. FRACT. 
(S) (KW/M2) (MJ/M2) (MJ/KG) (G/S) (KG/KG) (KG/KG) (KG/KG) 


186 138 0.45 0.00 0.041 620.8 0.14 0.000 
189 i 0.45 0.10 0.045 563.0 0.13 0.000 
a2 2.8 0.46 Ovs OL04952 51971 OFk2 0.000 
195 550 0.47 0.26 0.054 472.8 OFTl 0.000 
198 3.9 0.48 0782 0.057 454.9 0.10 0.000 
201 8 0.49 C787 0.059 434.9 0.10 0.000 
204 ag 0.50 0.40 0.061 419.4 0.11 0.000 
207 0.8 0.50 0.42 0.063 404.4 0.10 0.000 
210 576 OF51 0.44 0. 065¢% 3929S 0.10 0.000 
213 2.8 0.52 0.48 0.067 383.8 0.11 0.000 
216 Leo ls 0755 02,069.85 37058 0.10 0.000 
219 2.6 OF53 0.69 O.0%E Os 36292 0.10 0.000 
222 x 0.54 0796 O07, 07208 35557 0.09 0.000 
ae 4.4 0.56 136 0.074 346.0 0.09 0.000 
228 276 0.56 2.04 0.076 338.2 0.09 0.000 
231 6.2 0.58 3.07 0.07846 3325 0.08 0.000 
234 Shed 0.59 4.53 0. O78" 32785 0.08 0.000 
237 3.8 0.60 6.47 0.079.086 321-54 03077 0.000 
240 3.4 O362),@ 18.69 0.080 307.6 0.07 0.000 
243 8.0 0.64 1ia72% 0.08248 2902 0.09 0.004 
246 iGi.0 0.82 14.83 0.083 284.0 0833 0.019 
249 166.6 L732 18.02 0.084 276.7 0.72 0.040 
ha Vd 264.5 2d 21.08 0.087 266.1 0.99 0.046 
252 320.2 3.08 23279 0.0908 25835 ah 0.055 
258 344.9 4.11 26200 % 0.092" 25209 ree 0.059 
261 336.2 Sel 27.63 0.095 244.5 1.28 0.051 
264 nal 6.05 28368 % .0.097-6% 237226 1526 0.046 
267 301.6 6.96 29.520 = 0.099%8E 235471 1426 0.042 
270 288.8 7 382 293315 O07 TOWE 22652 ino 0.044 
273 2957.0 8 SAL 29X13; = 0: T0509 22404 Ls25 0.043 
276 29553 9.60 28 479.4 0.104. "5 21859 1.24 0.042 
219 PAM LA LOR49.D 28 37. 0.1068 218722 1.24 0.043 
282 301.0 1139.0 627 S78 OPO 21755 1.24 0.041 
285 29156 1227.8 27563.6 0, LO8e8 213.02 Lang 0.035 
288 27422 1309.8 27.38.% OSL08'92 210.5 Lets 0.035 
291 27672 L3U920 427 124.~ 051084 2152 DSi 0.040 
294 294.0 L480.8 27 5202 021094" 219 1.18 0.041 
27 296.6 LSMCOM 27 30." OPLT0S2 2076 Leo 0.036 
300 289.0 16%56.8 27 0.8 OSL10*"* 207.5 Lor 0.034 
303 284.7 LHGi1.8 275615 O04110%4 203-8 1. 0.037 
306 PAST IRT 18:80.8 2852h8 OFLIOVs 204.5 Lrg 0.038 
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309 
312 
ai) 
318 
S2r 
324 
B27 
330 
333 
336 
339 
342 
345 
348 
351 
354 
She) 
360 
363 
366 
369 
See 
5/5 
378 
381 
384 
387 
390 
oo3 
396 
599 
402 
405 
408 
411 
414 
417 
420 
423 
426 
429 
432 
435 
438 
441 
444 
447 
450 
453 
456 
459 
462 
465 


312% 
324. 
331. 
340. 
338. 
Sele 
Es 
343. 
so} 
348. 
368. 
B79. 
374. 
2} Oi ie 
Soe 
B50), 
5507, 
5565 
Sh a)e) 
368. 
364. 
348. 
342. 
345. 
50. 
347. 


363 


SM ye 
361. 
338. 
309. 
2935 
280. 
PAE 
247. 
253% 
216. 
209; 
198. 
76 
Lol 
Syke 
140. 
137 
141. 
AW oy. 
147. 
157. 
134. 
Las 
138. 
LESH 
L328 


FNYRMPUNDUDPWODONKYPNHMNWHYYNHDWYNOKFESFOWOWO AN FONAWAUFUONMFKFUWODWDAMAODOWOWOUOUNFEAHAOF AN 


Se) Soe ee Eero e @ jo) ©) ee) ) te) mele) ©) @ fe) eo) @) ©) ©) &) ee ee eyo) oo) () fe) Cas |) GIS ] © |] SS] &) 


. 109 
.109 
geiy) 
.110 
e109 
. 108 
LOT, 
LO, 
7107 
.105 
.104 
a.03 
.102 
- bor 
. 100 
.098 
.096 
:093 
O94 
.090 
.089 
.086 
.084 
.080 
.077 
: O75 
.072 
.069 
.068 
.065 
.063 
.060 
.057 
.054 
.O51 
.049 
.046 
044 
042 
.040 
.037 
.034 
. 032 
.029 
.027 
.025 
.022 
.020 
.018 
Oi? 
025 
.013 
OL 


205. 
205. 
206. 
204. 
201. 
204. 
208. 
207. 
206. 
2098 
2iif 
214. 
218% 
Zi 
220; 
229% 
228. 
230% 
241. 
244. 
250. 
254. 
263% 
2098 
282. 
2938 
305. 
32 
328% 
343. 
353% 
37 1g 
400. 
424. 
448. 
469. 
501m 
528. 
393 
584. 
6328 
684. 
741. 
808. 
873: 
952% 


OF RPHMYNONUODKFOFLFOAWNODOUNHAWNHOOKFNAWAWHKFODEFSFUWKYPNNOHKFPNIHDAwOOH 


1072. 
11973 
127% 
13986 
1549. 
1794. 
2098. 
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mNwWwWnN Fw NYC 


WMP FWWWWHDYNDNHNYNNNNNNMNNNNHNNNNNRPRP RRP PRP RP RP RP RP RP RPP RP RRP RP RP RP RP RP RP Pe 


(S) (@) (2) @)e) (2) te) () () (=) (&) ) fe) @) (2! (©) ) (ee) (2) (ej) (e) (ove) (e) (e) 2) (@) (=) fe) (@) ee )fe) (eo) (=) © So) (=) Ser) SI] IM &] OS] Oo Oo @ | 


.040 
.039 
.042 
.043 
.039 
.039 
041 
041 
.040 
.046 
.047 
045 
042 
.040 
042 
.040 
039 
.041 
.047 
.046 
045 
.047 
.047 
054 
.052 
O57, 
.060 
.061 
.060 
.057 
Rise ik 
.046 
.046 
.038 
.037 
.033 
034 
.035 
.030 
.025 
.023 
+028 
.022 
.026 
£029 
035 
41032 
.039 
041 
.040 
.042 
.043 
.043 


468 
471 
474 
477 
480 
483 
486 
489 
492 
495 
498 
501 
504 
507 
510 
513 
516 
54g 
522 
S20 
528 
yout 
534 
537 
540 


WWUWMNan oMn 0 


WMonomnwormonn-s 


BPWRONOOMORONPKYU:- 


SrOvorOrO O11 FOO (OOOO OOO (Ott OrOnO Oro ae 


(oe) (oe) (eel (ete) fe 2) fe) eM) ee 2) ee) OS eS) SSB] = 


.010 
.009 
.008 
.006 


004 


.003 
.002 
.001 
.001 
.001 
.001 
.001 
.001 
.001 
/O0F 
.001 
.001 
.001 
.001 
.001 
.001 
.001 
.001 
.001 
.001 
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(eWie) (e)(e) ee) (=) erate) eet ee ey) I OSS} 


(ej le} (ete) (e) [e) (eo) ey (@) (eee) (eo) (ej (oe) (a) (e) (= 


leV(eme Koka eo) (omel(o\(eo=) (oes (or (=) ee) (See Te Se) (=) 


047 
.045 
#035 
.036 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 
.000 


APPENDIX B: 


EXECUTING FURNITURE FIRE MODULES FOR TIME 


FOR ITER = 0 THE ICOUNT IS 3 
FOR ITER = 1 THE ICOUNT IS 2 
FOR ITER = 2 THE ICOUNT IS 1 


RESULT OF RADIATIVE HEAT TRANSFER ANALYSIS 


THE VAPOR PRESSURE FOR CO2 FOR LOWER AND UPPER GASES AND FOR H20 FOR LOWER AND 


UPPER GASES ARE1.428180332961E-6, 


5 .6448896955519E-3, 
FLAME.AREA RAD.FRACTION 


2*0. 


FLAME . HEIGHT 


1 0.4196E+01-0.3151E+00 0.1335E+01 0.3170E+00 
ELEM . AREA 


ELE SCAT/COSAV 
NETPOWER (KW) 

1 0.0000E+00 
0.3008E+03 

2 0.0000E+00 
0.8876E+02 

3 0.0000E+00 
0.2992E+00 


1 0.1000E+01 
0.2885E+01 
2 0.1000E+01 
0.2907E+01 
3 0.1000E+01 
0.2888E+01 
4 0.1000E+01 
0.2876E+01 
5 0.1000E+01 
0.2866E+01 
6 0.1000E+01 
0.2875E+01 
7 0.1000E+01 
0.2861E+01 
8 0.1000E+01 
0.2840E+01 
9 0.1000E+01 
0.2878E+01 
10 0.1000E+01 
0.3099E+01 
11 0.1000E+01 
0.3109E+01 
12 0.1000E+01 
0.3113E+01 
13 0.1000E+01 
0.3082E+01 
14 0.1000E+01 
0.3138E+01 


0 


0. 


.4196E+01 


6311E+03 


.6954E+03 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4840E-01 


.4840E-01 


.4840E-01 


.4840E-01 


.4840E-01 


EMISSIVITY BLCKBDYPOW 


0 


0. 


.59412E+00 0.1384E+03 


1920E+00 


.9883E-03 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


0. 


0 


1669E+01 


.4598E+00 


. 1361E+02 


. 1361E+02 


. 1361E+02 


.1361E+02 


. 1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


. 1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 
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250.000 SEC 


FLAME . RADIUS 


IRRADIANCE 


QO. 


0. 


0. 


5967E+01 


9359E+00 


8951E+00 


. 7810E+02 


. 7859E+02 


. 7818E+02 


. 7790E+02 


. 7769E+02 


.7789E+02 


.7756E+02 


.7710E+02 


.7794E+02 


.7764E+02 


. 7785E+02 


.7792E+02 


. 7730E+02 


. 7844E+02 


SAMPLE TABULAR OUTPUT OF FAST/FFM AT TIME 250 SECONDS 


RADIANCE 


0. 


0. 


1384E+03- 


1669E+01- 


.4598E+00 


.1361E+02 


.1361E+02 


. 1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


. 1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


15 0.1000E+01 
0.3166E+01 

16 0.1000E+01 
0.3056E+01 

17 0.1000E+01 
0.3082E+01 

18 0.1000E+01 
0.3147E+01 

19 0.1000E+01 
0.2847E+01 

20 0.1000E+01 
0.2857E+01 

21 0.1000E+01 
0.2906E+01 

22 0.1000E+01 
0.2925E+01 

23 0.1000E+01 
0.2900E+01 

24 0.1000E+01 
0.2932E+01 

25 0.1000E+01 
0.2906E+01 

26 0.1000E+01 
0.2864E+01 

27 0.1000E+01 
0.2905E+01 

28 0.1000E+01 
0.3185E+01 

29 0.1000E+01 
0.3108E+01 

30 0.1000E+01 
0.3094E+01 

31 0.1000E+01 
0.3174E+01 

32 0.1000E+01 
0.3134E+01 

33 0.1000E+01 
0.3083E+01 

34 0.1000E+01 
0.3146E+01 

35 0.1000E+01 
0.3081E+01 

36 0.1000E+01 
0.3056E+01 

37 0.8648E+00 
0.8761E+02 

38 0.9986E+00 
0.2091E+01 

39 0.9940E+00 
0.4194E+01 

40 0.9982E+00 
0.2139E+01 


.4840E-01 


-4840E-01 


-4840E-01 


.4840E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


.4473E-01 


-4473E-01 


.4840E-01 


.4840E-01 


.4840E-01 


.4840E-01 


.4840E-01 


.4840E-01 


.4840E-01 


.4840E-01 


-4840E-01 


. 3000E+03 


.6562E+01 


.8749E+01 


.6562E+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. LOOOE+01 


. 9000E+00 


. 9OO00E+00 


. 9000E+00 


. 9000E+00 


0. 


0. 


0. 


.1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


. 1361E+02 


.1361E+02 


. 1361E+02 


.1361E+02 


.1361E+02 


. 1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


. 1361E+02 


.1361E+02 


.1361E+02 


.1361E+02 


. 1361E+02 


.1361E+02 


.1361E+02 


5690E+00 


.5356E+00 


5356E+00 


5356E+00 
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.7902E+02 


.7675E+02 


.7729E+02 


. 7864E+02 


.7726E+02 


. 7749E+02 


.7858E+02 


.7901E+02 


. 7843E+02 


.7915E+02 


.7858E+02 


.7764E+02 


. 7854E+02 


.7942E+02 


.7783E+02 


. 7754E+02 


. 7920E+02 


. 7837E+02 


. 7730E+02 


. 7862E+02 


.7727E+02 


. 7675E+02 


.8935E+00 


.8896E+00 


. 1068E+01 


.8978E+00 


.1361E+02 
. 1361E+02 
. 1361E+02 
.1361E+02 
.1361E+02 
. 1361E+02 
.1361E+02 
.1361E+02 
~L361E+02 
.1361E+02 
.1361E+02 
. 1361E+02 
.1361E+02 
. 1361E+02 
.1361E+02 
.1361E+02 
.1361E+02 
.1361E+02 
.1361E+02 
.1361E+02 
.1361E+02 
.1361E+02 
.6014E+00 
.5710E+00 
. 5889E+00 


.5718E+00 


41 0.9968E+00 0. 
0.3586E+01 

42 0.9986E+00 0. 
0.3552E+01 

43 0.9986E+00 0. 
0.8773E+01 

44 0.9986E+00 0. 
0.3552E+01 

45 0.9975E+00 0. 
0.6000E+01 

46 0.7700E+00 0. 
0.1157E+03 


POWER IN KW ABSORBED BY ALL 


Time = 


Upper temp(K) 
Lower temp(K) 
Upper vol (m**3) 
Layer depth(m) 
Ceiling temp(K) 
Up wall temp(K) 
Low wall temp(K) 
Floor temp(K) 


Plume flow(kg/s) 
Pyrol rate(kg/s) 
Fire size(W) 


Plume in ul(W) 
Plume in 11(W) 
Vent fire(W) 

On target(W/m*2) 
Pressure (Pa) 


8749E+01 0.9000E+00 


2044E+02 0.9000E+00 
2725E+02 0.9000E+00 
2044E+02 0.9000E+00 
2725E+02 0.9000E+00 
3000E+03 0.9400E+00 


SURFACES 


250.0 seconds. 


412. 
300. 
144. 

0. 
318. 
els. 
304. 
300. 


Oo WM DW CO WM Oh Oo 


2.373E+00 
3.491E-02 
6.527E+05 
0 .000E+00 
0.000E+00 
6.527E+05 
0 .000E+00 
9.178E+00 
-1.644E-01 


Upper layer species 


N27 
O2ane zs 
co2 2a] 
CO ppm | 
TUHG) 20) 
H20 sez 
OD 1/m | 
Lower 
N22 | 
02 x | 
coz) Zig} 
CO ppm | 
TUHGs ez} 


1188 
19-6 
0.693 
SAGs 
ele) 
0.000 
2.46 


layer species 


7.903 

ZO 
7.898E-04 
0.432 
1.980E-03 


0 


0 


.5356E+00 


.4807E+00 


-4807E+00 


.4807E+00 


.4807E+00 


.4634E+00 


0 


0. 


.9910E+00 


.6738E+00 


.8384E+00 


.6738E+00 


. 7253E+00 


8737E+00 


1S345.1530768451 
FURNITURE BURNRATE IS 954.375025769 KW AT TIME 250. SEC 


300.0 


0 .000E+00 
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.5811E+00 


. 5|000E+00 


.5165E+00 


. |000E+00 


.5052E+00 


.4880E+00 


H20 % | 0.000 
OD 1/m | 3.708E-03 


FIRE SPREADING INFORMATION 
PANEL # IS1,# OF POLYGONS IS1 
PRINT POLYGON #, AREA, FLAME #, & FLAME FOOT 
be 024026 99Gley) b1 3E+2 
PRINT FIRE SPREAD VECTORS, MAGNITUDES, AND STATUS 
BaL7 , 108305 0191 S88969618srea0 ss, 

Al], © 9236954780 31915389696188a0y & 1 

TESS, IS69SRGTOT LOI SS89696IC ores 1 

Pear) 10-305, O. 191 5389696188888 0.5 1 

fois, 102305, "O19 1SS8969618 SREP 0-5 2 

PANEL # IS2,# OF POLYGONS IS1 

PRINT POLYGON #, AREA, FLAME #, & FLAME FOOT 

024356, 1, (076/00931109255 

PRINT FIRE SPREAD VECTORS, MAGNITUDES, AND STATUS 

7.83, 9.682975612144, 0.1989480023767, 

pel? 9.6829756121446 9% 001989660023767, 

7.17, 9.6254528219319980.28564365031173, 

Poea,  9.62945282198 182107 564565031173, 

7.83, 9.682975612144, 0.1989480023767, 

PANEL # 183, OF POLYGONS IS1 

PRINT POLYGON #, AREA, FLAME ##, & FLAME FOOT 

170-4026. (0.670162470819 

PRINT FIRE SPREA” VECTORS, MAGNITUDES, AND STATUS 

Perot ibe ozeso ee 669), 0219868786424833, 0., 3 

77197182824935, 10.305, 0.1988786424833) 0., 73 

7099660034720 587205305, 0.85636714632239) .0.,% 3 

7.09966003472550°95695; 028563671432239,..0., 3 

7215/182824935 e009 5695 ,@ 0.1988786424833;%.0., 83 

PANEL # 1IS4,# OF POLYGONS IS1 

PRINT POLYGON #, AREA, FLAME #, & FLAME FOOT 

1, 0.4356, 1, 0.6697461236551 

PRINT FIRE SPREAD VECTORS, MAGNITUDES, AND STATUS 

Pols, 1083030583052) 8 0. 199294989647TY0C0S, 

7 bo, 8 10, 92305830521 2 0. B99294989647 TR 2t0F, 

7.83, 10.37058109542, 0.8567834903877, 0O., 
Oy, 
Cr; 


ooo0oo 
NM BD fH 


7.17, 10.37058109542, 0.8567834903877, 
7.17, 10.31305830521, 0.1992949896471, 


FFE 


ELEMENT STATUS, TEMPERATURE (K), BNTIME(S), BURNFLUX(KW), & WALL FLAME 
DEPTH (M) 

1, 700., 38341.36999276, 651.4142349, 0. 

-1, 700., 65706.59054022, 268.7142386945, 0. 

1 700., 38415.30465874, 653.9898507037, 
700. 27847.75608523, 633.8904053881, 


Slo OF © 


a ) 

1, OO SAT BSE78 1007 73198-81649 [4995332359 ; 

tT, 700. 46827896 925919845 %-563427276318362 ; 

1, 700. 282229912069 76816720626 25538692056, 

1, 700., 22879.0818594, 624.0894625591, 0. 

1, 700 220G7154502473002),-85546:5030259223 0: 

1, 700., 16807.6914365, 512.2205860536, 7.6285999625314E-3 


127 


DP PRP PRP PRP RP RPP RPP PRP PRP PRP RP PRP RP RP RPP eB 


OONHA UNF WH FE 


700., 


.893863656785, 
901858700872, 
895982677837, 
-457429804298 , 

. 726190425738, 
.460030217096, 

. 394900311387, 
-404974470259, 
.440986126191, 

. 1378882971542, 
. 1382875864508, 
. 138378543403, 
. 3216153806969, 
. 325077916475, 
. 2334888332552, 
-4111304454612, 
4132264307185, 
4188314614041, 
21377632958955; 
. 2282681510027, 
. 139561409771, 
. 3277862224224, 
. 327281053455, 
. 3282383643718, 
4182379305449, 


0 
0 
0 
0 
0 
0 
0 
0. 
0 
0 
0 
0 
0 
0 
0 
0 


, 19375. 
819617. 
mL 16768. 
/ 15912. 
| 18394. 
mah 16234. 
, 16088. 
. 19888. 
, 27640. 
p37.96k. 
27976. 
27927. 

., 32098. 
WE? 7.938% 
, 16961. 
ee 74776 
, 23024. 
25360: 
, 15538. 
i ark63190 
20158. 
ep 22.9598 
ual 63274 
, 16574. 
6215. 

10659. 
30. 

, QCWALL 


468501, 

23141621, 
54628619, 
08532233, 
096372557 
21333122. 
10075563, 
63756089, 
97980879, 
95113684, 
S9L0Z1 23), 
72606525, 
86679194, 
10493834, 
19467181, 
68100926, 
41466151, 
44853421, 
87761288, 
38986986, 
872199725 
75568641, 
80847889, 
62612447, 
6620193, 
959153745 


, QRFLUX 


DO9 To) 210007 Zon 
57568758592109; 
511.744271116, 
529 .8676341817, 
545 .5033658545, 
509 .5772987183, 
520.9994631318, 
588 .8479224669, 
617 .9094057897, 
625 .9001048782, 
631.9459912086, 
635.1533990664, 
TOG 325992224795 
636.3141743875, 
517.3414304011, 
623 .9658431948, 
626 .4669165392, 
628 .4169400123, 
531. 160207338); 
517 .0406819988, 
629 .8621308906, 
625.0850027451, 
516 .5919270604, 
521.3818867327, 

518 57622530252) 
536 .4267594996, 


, SCALMS 


7.6090734896342E-3 


ae 


608324337309E-3 


1.3665331605362E-2 


NSRP RP PRP RPRPNYURDYN PRP RP Re 


.366708716641E-2 

.3658311258494E-2 
.7716516347934E-2 
.7717620314169E-2 
.7678453211916E-2 
.962297739792E-3 

.0680271618789E-3 
.9544533262175E-3 
.355945881842E-2 

.3439794416465E-2 
.3558959603248E-2 
.7717650514365E-2 
.7598350754424E-2 
.7686497724112E-2 
.5795685018843E-3 


7.6158223074965E-3 


Va 
ik 
1f 
Ll. 
iy: 


6133060780526E-3 
3641837683915E-2 
3629403794119E-2 
3656189604989E-2 
77089674442E-2 


1.7709574976736E-2 


a 


, SCALHT 


7603633878398E-2 


6.75 


8518425982 


64 .49031314678, 
98077620093, 
. 597043814858, 
. 28511526181, 
.07695533689, 
. 28085977568, 
94587482306, 
49187519249, 
.33115852878, 


NONNMNHNNHNHF WE 


.9200140983791E-4, 
.8094835364683E-5, 
.917029961399E-4, 

.2456458731578E-4, 
.0435195951577E-4, 
.2436500598888E-4, 
.3067543425569E-4, 
.3090630422256E-4, 
.2583337576681E-4, 


2 


6 
6 
6 
6 
6 
5 


. 663632239147 


6.779455950589 


.620780813278 
. 779521988194 
.629639978113 
. 571730388409 
576503485384 
. 70134730918 


64. 
64. 


63% 


63. 


65. 


65. 
64. 


64, 30953371725, 


64 .82964492049, 
.41033646555, 
. 13494432016, 
.67402831532, 
02559479421, 
.64552975078, 


64.97111689875, 


64.81990810553, 


02427864788, 
23854134671, 


6829372234, 


87565104909, 
39399459574, 
54102739649, 
96218069213, 
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2.3198881276931E-4, 
2.2993399881191E-4, 
2.2991778523806E-4, 
2 .3228932965644E-4, 
2.3196207717688E-4, 
2.3169066298695E-4, 
2.3208802406063E-4, 
2.3193849904426E-4, 
2.2954919183756E-4, 
2.2595848359531E-4, 
1.9228970711886E-4, 
2.2442145579029E-4, 
2.2465751967738E-4, 
2.1585954371332E-4, 
2.2457854809962E-4, 
2.3155331293433E-4, 


DNDNHNDMNMN WN 


5.439770585915 
6. 
6 .09728521708 

5 .443824335522 

5 .568484867426 

. 703487539841 
-447340254469 
. 5937407426096 
#16953,7033937. 
. 588693849189 
651158742362 
6 .644322198594 

6. 
7359749966737 
6. 
»: 


032964997961 


63541143549 


637311440089 
423979448524 


26, 
245 
28, 
“het 
30, 
oe 
Sai 
SG. 
34, 
aoe 
36, 


oOooOoooaooooffo°od 


5080200203823, 
4183382417381, 
.676763157653E-2, 
604456885221E-2, 
1377573564541, 
2338554252884, 
3250208026167, 
3216292475846, 
4184050004019, 
4131561803691, 
503926920194, 


64 .02580875025, 
64 .92970323311, 


63 .93088533261, 
65 .58800334798, 
64.76179633568, 
63 .69042700506, 
65.00711442756, 
632035/5982420501, 
63 .13780732048, 


65.80756989131, 
64.21977916703, 


2.254470181178E-4, 
2.3023948494752E-4, 


.3201011384592E-4, 
,3170326181462E-4, 
.3079576379553E-4, 
,3229090147856E-4, 
.3154709133589E-4, 
.3193980110939E-4, 
2.2509651021806E-4, 


NO NM NM NW fH fo 
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6 .600450349988 
6.570317615852 


2.2973849355817E-4, 6.555957260278 
2.3215697864253E-4, 5.63470605048 


.496419074609 
SOO 201407 29 
. 573834539398 
-494954026664 
463757387522 
. 514809114663 
5.728500906668 


WaAOnInnun 


APPENDIX C: SAMPLE OF PLOTTING FAST/FFM PREDICTIONS 


Explanations of generating these plots are in Section 7, Structure of FFM 
Output. 
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FURNITURE FIRE SPREAD 


Figure C-1. Three-dimensional plot of FFM flame spread data 
in Unit 13 generated by selecting IPLOT = 3 
(See Table 7-4) 


BURN AREA FRACTION 


Furniture Calorimeter 
0 seat 
A back rest 
x left arm rest 
oO right arm rest 


Y . haieed 


FAST/FFM 
— seat 
-- back, left, 

and right rests 


22 


WO eCeescerccas 
ee ec ere cence 
° 


Peccecce 
ee e2ee 


— 


eresesreszoccce 


@.0¢ x 
0.00 100.00 200.00 


TIME (SEC) 


Figure C-2. 


300.060 400.00 


Burn area comparisons to full scale data for 


the HO/FRPU furniture mockup. 
data obtained from Unit 12 file (see Table 7-3) 


FFM predicted 


BURNRATE (kW) 


1200.0 


Figure C-3. 


o Furniture calorimeter 
— FAST/FFM 


100.00 | 200.00 300.08 400.00 


TIME (SEC) 


FFM burn rate predictions compared to full 
scale data for the HO/FRPU furniture mockup. 
FFM predicted data from Unit 10 file 

(see Table 7-1) 


MASS LOSS RATE (g/s) 


40.02 


28 .8@ 


10. 0¢ 


Figure C-4. 


o Furniture calorimeter 
— FAST/FFM 


300.00 


400.80 


TIME (SEC) 


FFM fuel loss rate predictions compared to full 
seale data for the HO/FRPU furniture mockup. 
FFM predicted data obtained from Unit 10 file 
(see Table 7-1) 
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SOOT EXTINCTION AREA (m2/kg) 


2088 . OC 


1500. 0¢ 


1000 .0@ 


Figure C-5. 


------ Furniture calorimeter 
FAST/FFM 


100.00 200.60 300.00 400.00 


TIME (SEC) 


FFM soot extenction area predictions compared 
to full scale averaged value for the HO/FRPU 
furniture mockup. FFM predicted data obtained 
from Unit 11 file (see Table 7-2) 


LoD 


PRODUCT FRACTION (%) 


$00.02 


250. 02 


-250.0 


@.00 100.00 280.00 300.20 420.02 


TIME (SEC) 


Figure C-6. FFM predictions of furniture mockup products. 
Plotting data obtained fron Unit 10 file (see 
Table 7-1) 
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CONVECTIVE HEAT RATE (kW) 


4.02 


Figure C-7. 


| 
100.00 200.88 308.80 400.00 


TIME (SEC) 


FFM predictions of total convective heat rate 
on the furniture mockup surfaces. Plotting 
data obtained fron Unit 10 file (see Table 7-1) 
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ZONE RADIATION POWER (kW) 


600.0 


400.0 


a) 
®& 
® 
® 


e 


-200.@ 


@.00 


Figure C-8. 
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i 
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t 
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108.00 206.89 300.88 400.00 
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FFM predictions of net radiative power output 
of furniture fire. Plotting data obtained from 
Unit 11 file (see Table 7-2) 


138 


TEMPERATURE (°K) 


450.0 


400.0 


350.0 


302.2 


Figure C-9. 


——upper gas layer 
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upper sre Bia eee 
wal] - eoreceme a Sores 
ae - Dein W all eee 
: utt eer ah ae Sand 
Hiei: gas layer floor 
180.00 200.80 300.08 400.00 


TIME (SEC) 


FAST predictions of gas zones and room 
temperatures. Plotting data obtained from Unit 
11 file (see Table 7-2) 
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SCALED HEAT RELEASE RATE (-) 


4.0 


2.8 


8.08 


Figure C-10. 


28.80 49.80 60.90 


SCALED TIME (MJ/m°) 


Scaled heat release flux versus time for HO 
fabric/FRPU foam material analoous to Fig. 3- 
3b. Plotting data obtained frm Unit 14 file 


(see Table 7-5) 
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SCALED FUEL RELEASE RATE (g/kJ) 


8.1 


0.1¢ 


8.85 


0.00 20.60 40.00 60.00 


SCALED TIME (MJ/m¢) 


Figure C-11. Scaled fuel release flux versus time for HO 
fabric/FRPU foam material analogous to Fig 3- 
4b. Plotting data obtained from Unit 15 file 
(see Table 7-5) 
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